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UNIT  OF  MEASURE 

ABBREVIATIONS  USED  IN  THIS  REPORT 

ft 

foot 

mi/h          mile  per  hour 

h 

hour 

min          minute 
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kilogram 
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liter 

pet          percent 
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meter 

COMPUTER   SIMULATION   OF   OXYGEN    REQUIREMENTS    FOR    EVACUATING   MINERS 

By  John  C.  Edwards     and  Glenn  N.  Grannemann 


ABSTRACT 

The  Bureau  of  Mines  developed  a  mathematical  model  to  predict  the  oxy- 
gen requirements  for  a  miner  evacuating  a  mine.  The  model  accounts  for 
variations  in  the  miner's  escape  speed  with  mine  airway  height,  and  for 
the  empirically  determined  expressions  for  the  weight-specific  oxygen 
uptake  rate  for  upright  walking,  stooped  walking  (duck-walking) ,  and 
crawling.  A  Fortran  computer  program  that  is  user  oriented  was  prepared 
based  upon  the  model's  constitutive  equations.  Specific  instructions 
for  utilizing  the  program,  as  well  as  applications  to  miner  evacuation 
from  a  mine,  are  presented. 

INTRODUCTION 

Federal  regulations  require  that  every  person  who-  enters  an  under- 
ground coal  mine  in  the  United  State  be  supplied  with  a  self-contained 
self -rescuer  (SCSR).  An  SCSR  is  required  to  supply  oxygen  to  the  user 
for  at  least  60  min.  In  support  of  this  requirement,  the  Bureau  of 
Mines  has  directed  research  into,  and  provided  evaluation  of,  SCSR  tech- 
nology (1-2). ^  To  evaluate  the  most  efficient  utilization  of  an  SCSR  in 
miner  withdrawal  and  rescue,  it  is  necessary  to  have  an  adequate  model 
of  the  miner's  demand  on  oxygen  while  traveling  through  mine  airways. 
For  this  purpose  the  Bureau  developed  a  Fortran  computer  program  that 
evaluates  the  oxygen  consumption  at  speed-dependent  oxygen  uptake  rates 
for  miners  traveling  over  specified  paths  in  the  mine.  This  model  uti- 
lizes a  generalization  of  the  dependence  of  the  weight-specific  oxy- 
gen uptake  rate,  Vq2»  upon  the  miner's  egress  speed  and  posture.  This 
forms  the  first  constitutive  equation  of  the  model. 

The  miner's  average  speed  in  a  mine  airway  has  been  shown  (3)  to  de- 
pend upon  the  airway  height.  From  Berry's  data  (_3)  a  second  constitu- 
tive equation  was  formed  for  the  model.  The  two  constitutive  equations 
form  the  core  of  the  Fortran  program  that  evaluates  oxygen  consump- 
tion by  a  miner.  The  equations  selected  for  this  model  are  discussed 
in  the  following  section.  In  the  section  on  Data  Preparation,  the  pro- 
gram utilization  is  described.  In  the  section  on  Application  the  model 
is  applied  to  the  evacuation  of  32  miners  from  a  mine  that  contains 
41  airways. 


1  Physicist,  Pittsburgh  Research  Center,  Bureau  of  Mines,  Pittsburgh,  PA. 
2Systems  analyst,  Product  Research  Inc. 

3Underlined  numbers  in  parentheses  refer  to  items  in  the  list  of  references  preced- 
ing the  appendixes. 


MODEL 


The  weight-specific  oxygen  uptake  has 
been  determined  to  be  statistically  best 
represented  (4-7)  by 


V0  =  7.06  +  0.001  S2 


(1) 


for  walking  in  an  upright  position  and  by 


v02  =  4*5  +  °*2  S 


(2) 


for  running, 


ro2 


Vn.  =  oxygen  uptake  rate,  mL  02 
per  kilogram  body  weight 
per  minute 


and 


S  =  miner's  travel  speed, 
m/min. 


Either  the  parameters  aj ,  0j,  Xjj  for  i  = 
1,2  and  j  =1,2  are  supplied  by  the  pro- 
gram user,  or  the  program  uses  values 
from  equations  1  and  2  as  default  values. 
For  airway  heights  significantly  less 
than  the  miner  height,  the  miner  will 
have  to  either  duck -walk  or  crawl.  In 
either  case  the  oxygen  uptake  will  in- 
crease if  the  miner  attempts  to  maintain 
his  walking  speed.  Based  upon  Berry's 
data  O,  table  15)  ,  the  following  speed- 
dependent  oxygen  uptakes  can  be  postu- 
lated for  crawling  and  duck  walking, 
assuming  a  linear  dependence  upon  speed: 

Vq2  =  0.81  S   (crawling)       (6a) 

V0  =  0.66  S   (duck-walking)   (6b) 


Equation  2  is  valid  for  speeds  in  excess 
of  134  m/min  (5  mi/h).  These  expressions 
have  been  utilized  in  studies  by  Kamon 
(_7-8) .  An  expression  similar  to  equation 
2,  but  more  general,  was  developed  by 
Taylor  (9)  for  many  species  of  animals  in 
locomotion.  Taylor's  expression  is  given 
by 


These  expressions  are  generalized  to  in- 
clude weight  dependence,  as  was  done  for 
walking  and  running  in  equations  5a  and 
5b: 

X3 
V02  =  a3M  S   (crawling)       (7a) 

Vq2  =  ot4M  4S   (duck-walking)   (7b) 


V02  =  18.0M-0.303  +  o.533M-°.316  S,   (3) 

where   M  =  the  animal's  mass,  kg. 

For  an  87-kg  miner   (an  average  50th- 
percentile  miner) ,  equation  3  reduces  to 


V0  =  4.65  +  0.13  S, 


(4) 


which  is  in  close  agreement  with  equation 
2.  To  increase  the  scope  of  the  computer 
program,  the  first  constitutive  equation 
of  the  model  was  written  for  upright  lo- 
comotion as 


Based  upon  Berry's  data  (3),  duck- 
walking  is  assumed  to  occur  wherever  the 
airway  height  is  less  than  50  in,  and 
crawling  wherever  the  airway  height  is 
less  than  30  in. 

Equations  5  and  7  form  the  first  con- 
stitutive equation. 

The  second  constitutive  equation  for 
the  model  quantifies  the  dependence  of 
the  average  miner  speed,  S,  on  the  airway 
height,  H.  Using  Berry's  data  (3,  fig. 
ES-1)  as  a  model,  the  following  constitu- 
tive equation  was  constructed  using  a 
linear  regression  computer  program: 


V0     =  aiM       S2   +  6,M 


12 


(5a) 


S  =   0.3048   (A+BH+CH2+DH3) 


(8) 


where   S<134  m/min  (walking),  or 


A2 1       ' 
a2M   S  +  6?M 


22 


v02  ~  a2"    o  -r  p2i 
where   S>134  m/min  (running). 


The  coefficients  A,  B,  C,  and  D  were 
determined  from  the  regression  analy- 
(5b)  sis  for  a  miner  walking,  crawling,  or 
duck-walking.  The  entry  height  H  is  in 
inches.   Airway  heights  and  lengths  have 


units  of  inches  and  feet  to  correspond 
more  readily  to  available  mine  data, 
whereas  physiological  data  are  in  metric 
units,  which  are  more  common  in  the  lit- 
erature. Values  of  the  coefficients  A, 
B,  C,  and  D  are  shown  in  table  1.  The 
program  user  may  either  supply  a  new  set 
of  coefficients,  or  the  program  will  use 
the  values  in  table  1  as  default  values. 

TABLE  1.  -  Escape  speed  equation 
coefficients 

Values 

Coefficient: 

A 252.2 

B -16.14 

C 408 

D -.00244 

The  mine  network  is  described  by  air- 
ways of  specified  length  and  height  that 
intersect  at  junctions.  Each  miner  is 
initially  at  a  prescribed  junction.  The 
miner  progresses  along  a  prescribed  path 
consisting  of  airways  until  the  "escape 
junction"  (the  mine  exit)  is  reached. 
The  program  also  permits  rest  periods  at 
specified  junctions.  The  miner  escape 
speed  along  each  airway  is  determined  by 
equation  8,  and  the  weight-specific  oxy- 
gen uptake  is  determined  from  the  escape 
speed  using  equations  5-7.  The  total 
oxygen  requirements  of  the  miner  are  de- 
termined as  a  function  of  elapsed  time 
and  distance  traveled. 

It  can  be  shown  that  the  volume  V 
(liters)  of  oxygen  consumed  by  a  miner  of 
mass  M  (kilograms)  after  traveling  for  a 
time  T  (minutes)  is  given  for  a  miner 
walking  by 


V  = 


M 


1,000 


cciM 


12, 


+  ptM   T 


(9a) 


and  for  a  miner  running  by 
M 


V  = 


1,000 


Xo  i        A22 
a2M   X  +  62M   T 


(9b) 


where  X  is  the  total  distance  traversed, 
in  meters. 


Equation  7b  is  readily  evaluated  for  a 
single  miner  running.  However,  for  a 
miner  traversing  a  path  at  variable 
speeds,  and  for  a  number  of  miners  tra- 
versing different  escape  paths,  a  com- 
puter program  has  a  definite  advantage  in 
the  evaluation  of  total  oxygen  demand 
made  upon  an  SCSR.  Aside  from  the  compu- 
ter program's  calculated  oxygen  require- 
ment for  a  miner  egressing  from  a  mine,  a 
compressed-oxygen  SCSR's  oxygen  supply 
will  never  be  less  than  the  constant  sup- 
ply rate  of  the  SCSR.  This  must  also  be 
considered  by  the  program  user.  The 
response  of  an  SCSR  to  a  breathing  demand 
that  varies  as  a  sinusoidal  function  of 
time  has  been  treated  elsewhere  (10)  as  a 
separate  mathematical  model. 

For  the  case  of  a  miner  walking  at  a 
constant  speed,  equation  9a  can  be  evalu- 
ated as 


V  = 


M 


1,000 


Xi  1    „ 
o,M        X2/T 


61M     ZT 


(10) 


For  a  constant  quantity  of  oxygen  avail- 
able, equation  9b  shows  a  linear  rela- 
tionship between  distance  traveled  and 
elapsed  time,  whereas  equation  10  shows  a 
nonlinear  relationship. 

For  a  specified  escape  distance,  it  is 
possible  to  evaluate  from  equations  9b 
and  10  the  influence  of  travel  speed  upon 
oxygen  requirements.  For  example,  an  87- 
kg  miner  will  walk  2,000  m  at  a  speed  of 
67  m/min  (2.5  mi/h)  in  29.8  min,  and  can 
run  the  same  distance  at  a  speed  of  147 
m/min  (5.5  mi/h)  in  13.6  min.  The  miner 
will  require  29.9  L  02  in  the  walking 
mode  and  40.1  L  02  in  the  running  mode. 
This  is  an  important  consideration  for 
planning  miner  rescue  with  an  SCSR  be- 
cause the  total  oxygen  requirement  is  as 
important  a  consideration  as  minimizing 
the  escape  time. 

There  are  two  limitations   imposed  upon 

a  miner's  breathing.   One  is  the  absolute 

(maximum  weight-specific)   oxygen  uptake, 

Va  „,„,   which  is   characteristic  of   a 
u2  ma  x ' 

miner,  and  the  other  is  the  duration  a 
miner  can  function  at  a  working  capacity, 
which  is   a  specified  fraction  of   the 


individual's  aerobic  capacity.  Bonjer 
(11)  has  developed  an  analytic  expres- 
sion that  relates  an  individual's  working 
capacity,  or  oxygen  uptake,  to  the  indi- 
vidual's aerobic  capacity,  or  absolute 
oxygen  uptake  over  the  maximum  endurance 
time  t  (min) : 


log  (T)  =  3.76  -  3(V0/V0,max) 


(ID 


For   this  model,   the   oxygen  uptake 
in   equation   11   is   replaced   by   a 


time-weighted  average  of  the  oxygen  up- 
take specific  to  each  airway  through 
which  the  miner  traverses.  The  program 
does  not  permit  the  miner  to  exceed 
either  the  constraint  of  an  absolute  oxy- 
gen uptake  or  an  endurance  time  defined 
by  equation  11.  Equation  11  shows  that  a 
miner  can  be  expected  to  endure  at  50  pet 
of  aerobic  capacity  for  ~3  h,  and  at  65 
pet  of  aerobic  capacity  for  ~1  h. 


DATA  PREPARATION 


In  preparation  for  program  usage,  an 
input  data  file  must  be  prepared  that  de- 
scribes the  mine,  the  miners,  and  the  es- 
cape paths.  Each  file  is  created  as  an 
ASCII  file  that  is  translated  by  a  pro- 
gram called  TRANS. FOR,  listed  in  appendix 
B,  into  a  binary  data  file.  The  binary 
data  files  are  read  by  the  Fortran  pro- 
gram RESCUE. FOR,  listed  in  appendix  C, 
which  performs  the  model  calculations. 

The  requirements  of  the  input  data 
files  are  described  under  the  headings 
that  follow.  Each  element  in  the  data 
file  is  written  in  free  format. 

A.  MINE  DESCRIPTION 

The  input  data  for  a  description  of  the 
mine  is  entered  into  the  ASCII  data  file 
XXXAIR.ASC,  where  XXX  is  a  prefix  as- 
signed by  the  program  user  to  distinguish 
mine  descriptions.  The  same  prefix  is 
used  for  the  files  described  below  for 
the  miners  and  the  escape  paths.  The 
data  for  the  number  of  airways  (NA) 
is  entered  into  NA  rows  of  the  file 
XXXAIR.ASC.  Each  row  must  have  the  re- 
quired elements  shown  in  table  2. 

TABLE  2.  -  Mine  airway  data  file 

Airway  Junction  Junction  Height,  Length, 
1 2 iri ft 

Junctions  1  and  2  designate  the  junctions 
that  terminate  the  airway.  The  data  file 
is  terminated  by  airway  999. 


B.  MINER  DESCRIPTION 

The  input  data  for  a  description  of  the 
number  of  miners  (NM)  is  described  by  the 
ASCII  data  file  XXXMIN.ASC.  There  are 
2  NM  rows  of  data,  where  two  adjacent 
rows  are  used  to  describe  each  miner. 
The  first  row  of  each  pair,  as  shown  in 
table  3,  lists  the  miner  identification 
number,  the  junction  from  which  the  miner 
starts,  the  identification  number  of  the 
escape  path,  and  the  absolute  oxygen  up- 
take of  each  miner.  The  second  row  of 
each  pair  permits  the  specification  of  as 
many  as  three  rest  locations,  and  the 
corresponding  rest  time.  The  end  of  the 
data  is  signified  by  the  assignment  of 
miner  number  999. 

TABLE  3.  -  Miner  data  file 

Miner  Start   Escape   Absolute  oxygen 
junction  path  uptake,  mL/(kg*min) 

C.  ESCAPE  PATH  DESCRIPTION 

The  input  data  file,  XXXPAT.ASC,  de- 
scribes the  escape  paths  through  an  enu- 
meration of  the  constitutive  airways  of 
each  path.  Each  path  can  contain  as  many 
as  seven  airways,  as  shown  in  table  4. 
If  additional  airways  are  required,  then 
a  new  path  number  is  specified  in  the 
column  titled  "Next  path."  Wherever  no 
additional  airways  or  paths  are  required, 
a  zero  is  entered.  To  mark  the  end  of 
the  data  file,  path  999  is  entered. 


TABLE  4.  -  Escape  path  data  file 


Path   Airway 
1 


Airway   Next  path 
7 


There  are  several  questions  which  the 
programs,  TRANS  and  RESCUE,  ask  the  user 
in  the  interactive  mode.   The  protocol 


for  these  questions  is  shown  in  appen- 
dixes B  and  C.  These  questions  relate  to 
the  selection  of  an  output  printing  de- 
vice, whether  or  not  rest  junctions  are 
to  be  included,  and  the  constants  for 
equations  5-7.  The  user  has  the  choice 
to  enter  new  values  or  use  the  program's 
default  values. 


APPLICATION 


The  program  was  applied  to  a  mine  with 
41  airways  and  26  junctions,  as  shown  in 
figure  1.   Junction  26  represents  the 


mine  portal  to  be  used  for  escape. 
Thirty-two  miners  were  located  throughout 
the  network. 
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FIGURE  1.  -   Schematic  of  mine  network. 
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The  three  ASCII  files  that  describe  the 
mine,  the  miners,  and  the  escape  paths 
are  listed  in  appendix  D.  A  sample  out- 
put for  miners  1  and  32  is  displayed  in 
appendix  E.  Upon  completion  of  travel 
along  an  airway,  the  airway  number  is 
listed  in  the  output,  as  well  as  the 
elapsed  time  and  distance  traversed. 
Additional  information  is  provided  re- 
garded the  miner's  cumulative  oxygen  con- 
sumption, weight-specific  oxygen  uptake, 
speed,  and  maximum  oxygen  uptake.  After 
the  miner  reaches  the  exit  portal,  the 
maximum  endurance  time,  t,  as  derived 
from  equation  11,  is  listed. 

In  appendix  D  the  computer  program  out- 
put for  miners  1  and  32  is  presented. 
These  two  miners  have 
78.6  kg,  and  absolute 
34.10  mL/(kg»min),  and 
distance,  3,000  ft,  to 
separate  paths  from 


identical  masses, 
oxygen  uptakes , 
travel  the  same 
the  exit  along 

distinct  initial 


junctions.  The  paths  differ  in  the  air- 
way height  along  various  segments,  there- 
by requiring  different  modes  of  locomo- 
tion, namely  crawling,  duck-walking,  and 
walking.  Miner  1  starts  at  junction  1 
and  proceeds  along  airways  2,  4,  5,  and 
41.  Miner  32  starts  at  junction  25  and 
proceeds  along  airways  40,  38,  35,  25, 
41.  Miner  1  reaches  the  exit  in  31  min, 
and  miner  32  reaches  the  exit  in  15  min. 
Total  oxygen  requirements  are  39.9  L  for 
miner  1  and  26.4  L  for  miner  32.  Average 
oxygen  costs  per  kilogram  mass  per  meter 
traveled  are  0.55  mL  02  for  miner  1  and 
0.37  mL  02  for  miner  32.  This  illus- 
trates that,  because  of  the  various  re- 
quired modes  of  travel,  a  miner  requires 
a  quantity  of  oxygen  during  the  escape 
that  is  dependent  not  only  upon  the  total 
distance  traveled,  but  also  upon  the  mode 
of  travel;  i.e.  ,  posture. 


CONCLUSION 


A  user-oriented  Fortran  computer  pro- 
gram was  developed  that  enables  those 
charged  with  planning  mine  safety  to 
rapidly  evaluate  the  oxygen  requirements 
of  miners  in  the  evacuation  of  a  mine. 
The  program  accounts  for  the  variation  of 
miner  egress  speed  with  airway  height, 


and  the  dependence  of  weight-specific 
oxygen  uptake  upon  miner  speed  and  pos- 
ture. The  program  can  be  used  in  an  ac- 
tual mining  situation  to  estimate  the 
time  for  miner  evacuation,  and  the  ade- 
quacy of  the  oxygen  supply  from  an  SCSR. 
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APPENDIX  A. — LIST  OF  SYMBOLS 

H  airway  height,  in 

M  miner's  mass,  kg 

S  miner's  travel  speed,  m/min 

T  miner's  travel  time,  min 

V  oxygen  volume,  L 

V0  oxygen  uptake,  mL/(kg»min) 

V0  max  absolute  oxygen  uptake,  mL/(kg»min) 

X  distance,  m 

Coefficients  in  oxygen  uptake: 

a] ,  61  walking 

a2,  82  running 

03  crawling 

04  duck-walking 

Exponents  of  miner's  mass  in  oxygen  uptake: 

X  1  t ,  A21  walking 

Aj2»  ^22  running 

X3  crawling 

\4  duck-walking 

t  maximum  endurance  time,  min 


APPENDIX  B.— LISTING  OF  TRANS. FOR 


SAMPLE  TERMINAL  DIALOGUE  FOR  PROGRAM  TRANS 
!  >>  COMMENTS  ADDED  LATER 
$  >>  SYSTEM  PROMPT 


$  RUN  TRANS 
ENTER  PREFIX  FOR  DATA  FILE  NAMES:   MNl 

PROCESSING  AIRWAY  FILE 

PROCESSING  MINER  FILE 

PROCESSING  PATH  FILE 

MINE  DATA  FILE  TRANSLATED 
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0001 

C 

0002 

C 

0003 

C 

0004 

C 

0005 

C 

0006 

C 

0007 

C 

0008 

C 

0009 

C 

0010 

C 

0011 

C 

0012 

c 

0013 

c 

0014 

c 

0015 

c 

0016 

c 

0017 

c 

0018 

c 

0019 

c 

0020 

c 

0021 

c 

0022 

c 

0023 

c 

0024 

c 

0025 

c 

0026 

c 

0027 

c 

0028 

c 

0029 

c 

0030 

c 

0031 

c 

0032 

c 

0033 

c 

0034 

c 

0035 

c 

0036 

c 

0037 

c 

0038 

c 

0039 

c 

0040 

c 

0041 

c 

0042 

c 

0043 

c 

0044 

c 

0045 

c 

0046 

c 

0047 

c 

0048 

c 

0049 

c 

0050 

c 

0051 

c 

0052 

c 

0053 

c 

0054 

c 

0055 

c 

0056 

c 

0057 

c 

PROGRAM  TO  TRANSLATE  ASCII  FILES  OF  MINE  RESCUE  DATA 

TO  BINARY  FILES.   THE  EDITOR  IS  USED  TO  CREATE  THE 

ASCII  VERSION  OF  THE  DATA  FILES.   THESE  FILES  ARE 

THEN  RUN  THROUGH  THIS  PROGRAM  AND  A  BINARY  VERSION  IS  WRITTEN. 

THERE  ARE  3  ASCII  FILES  USED  BY  THIS  PROGRAM. 

THESE  FILES  ARE  (THE  FIRST  3  CHARACTERS  OF  THE  NAME  INDICATE 
TO  WHICH  MINE  THE  FILE  REFERS.) 

MN1AIR.ASC 

THIS  IS  AN  AIRWAY  LIST  FOR  MINE  1.   IT  GETS  TRANSLATED  BY  THE 
PROGRAM  TRANS  INTO  A  BINARY  FILE.   THE  FORMAT  IS: 


AIRWAY   JUN      JUN 
NUMBER   1        2 


HEIGHT   LENGTH 
(IN.)    (FT.) 


ENTER  0  FOR  REMAINING  JUNCTIONS  IF  THERE  ARE  LESS  THAN  6 
JUNCTIONS  IN  AN  AIRWAY. 

AIRWAY  #  9  99  STOPS  THE  PROGRAM. 

THIS  IS  A  MINER  LIST  FOR  MINE  1.  IT  GETS  TRANSLATED  BY  THE 
PROGRAM  TRANS  INTO  A  BINARY  FILE.  THE  FORMAT  FOR  THE  FIRST 
LINE  IS: 


MINER 
NUMBER 


START 
JUNCTION 


ESCAPE 
PATH 


MAX  0-2 
INTAKE 


MINER  #  999  STOPS  THE  PROGRAM. 


THIS  IS  A  MINER  MASS  LIST  FOR  MINE  1.   IT  GETS  TRANSLATED 
BY  THE  PROGRAM  TRANS  INTO  A  BINARY  FILE.   THE  FORMAT  FOR 
THE  SECOND  LINE  IS: 


MINER    MINER 
NUMBER   MASS 


REST  REST 
JUNCT  TIME 
1        1 

MN1PAT.ASC 


REST  REST  REST 
JUNCT  TIME  TIME 
2        2        3 


THIS  IS  AN  ESCAPE  PATH  LIST  FOR  MINE  1.  IT  GETS  TRANSLATED 
BY  THE  PROGRAM  TRANS  INTO  A  BINARY  FILE.   THE  FORMAT  IS: 


PATH  AIRWAY  AIRWAY  AIRWAY  AIRWAY  AIRWAY  AIRWAY  AIRWAY  NEXT 
NUMBER  12      3        4       5       6       7     PATH 

PATH  #  999  STOPS  THE  PROGRAM. 

THESE  ASCII  FILES  ARE  TRANSLATED  INTO  A  DATA  BASE  FILE 
MN1DAT.DAT.   THE  .DAT  FILE  HAS  FOLLOWING  STRUCTURE: 

EACH  RECORD  IS  32  BYTES  (8  LONG  WORDS)  LONG. 
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0058 

C 

0059 

C 

0060 

C 

0061 

C 

0062 

C 

0063 

C 

0064 

C 

0065 

C 

0066 

C 

0067 

0068 

C 

0069 

C 

0070 

c 

0071 

0072 

0073 

0074 

0075 

0076 

c 

0077 

c 

0078 

c 

0079 

0080 

0081 

0082 

c 

0083 

c 

0084 

c 

0085 

0086 

c 

0087 

c 

0088 

c 

0089 

0090 

0091 

0092 

c 

0093 

c 

0094 

c 

0095 

0096 

0097 

0098 

0099 

c 

0100 

c 

0101 

c 

0102 

0103 

c 

0104 

c 

0105 

c 

0106 

0107 

c 

0108 

c 

0109 

c 

0110 

0111 

c 

0112 

0113 

c 

0114 

c 

50  RECORDS  FOR  LINE  1  DATA  FROM  MN1MIN.ASC 
50  RECORDS  FOR  LINE  2  DATA  FROM  MN1MIN.ASC 
100  RECORDS  FOR  DATA  FROM  MN1AIR.ASC 
100  RECORDS  FOR  DATA  FROM  MNlPAT.ASC 

PROGRAM  TRANS 

VARIABLES  FOR  FILE  NAMES 

BYTE  FILAIR(12) 
BYTE  FILDAT(12) 
BYTE  FILMIN(12) 
BYTE  FILNAM(12) 
BYTE  FILPAT(12) 

VARIABLES  FOR  ASCII  CHARACTERS 

BYTE  IDOT , ICHA, ICHC , ICHD , ICHI , ICHJ 
BYTE  ICHM , ICHN , ICHP , ICHR, ICHS , ICHT 
BYTE  ICHU 

TWO  BYTE  FILL  CHARACTER 

INTEGER  *2  IFILL2 

VARIABLES  TO  INDEX  AIRWAYS,  JUNCTIONS,  MINERS,  AND  PATHS 

INTEGER  *2  IAIR 
INTEGER  *2  IMIN 
INTEGER  *2  IPAT 

VARIABLES  FOR  LUNS 

INTEGER  *2  FAIR 

INTEGER  *2  FDAT 

INTEGER  *2  FMIN 

INTEGER  *2  FPAT 

FOUR  BYTE  INTEGER  FILL  CHARACTER 

INTEGER  *4  IFILL4 

FOUR  BYTE  REAL  FILL  CHARACTER 

REAL  *4  FILL 

INPUT  VARIABLES 

REAL  *4  HEIGHT, LENGTH  i AIRWAYS 

INTEGER  *2  NXTPAT  1NEXT  PATH  (IF  NECESSARY) 

SET  UP  ASCII  CHARACTERS  TO  USE  WITH  FILE  NAMES 
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0115 

C 

0116 

0117 

0118 

0119 

0120 

0121 

0122 

0123 

0124 

0125 

0126 

0127 

0128 

0129 

C 

0130 

C 

0131 

c 

0132 

0133 

0134 

0135 

c 

0136 

c 

0137 

c 

0138 

0139 

0140 

0141 

0142 

c 

0143 

c 

0144 

c 

0145 

1 

0146 

c 

0147 

c 

0148 

c 

0149 

0150 

0151 

0152 

0153 

100 

0154 

c 

0155 

c 

0156 

c 

0157 

c 

0158 

c 

0159 

0160 

0161 

0162 

0163 

0164 

0165 

0166 

0167 

c 

0168 

c 

0169 

c 

0170 

0171 

I  DOT 

= 

• 

ICHA 

= 

•A' 

ICHC 

= 

•c 

ICHD 

= 

'D' 

ICHI 

= 

•I' 

ICHJ 

= 

'J* 

ICHM 

= 

'M' 

ICHN 

= 

'N' 

ICHP 

= 

-p. 

ICHR 

= 

»R» 

ICHS 

= 

'S' 

ICHT 

= 

1  rp  1 

ICHU 

= 

'U' 

NULL  FILL  CHARACTERS 

IFILL2  =  0 
IFILL4  =  0 
FILL  =  0. 

SET  UP  LUNS 


FAIR 

=  7 

FDAT 

=  8 

FMIN 

=  9 

FPAT 

=  10 

ILUN  FOR  AIRWAY  FILE 

ILUN  FOR  BINARY  FILE 

!LUN  FOR  MINER  FILE 

ILUN  FOR  PATH  FILE 


GET  FILE  NAME  FOR  DATA  INPUT 

CALL  PRMPTA( '  ENTER  PREFIX  FOR  DATA  FILE  NAMES:   ',FILNAM,M) 

TRANSFER  PREFIX  OF  FILE  NAMES  INTO  REQUIRED  VARIABLES 

DO  100  I  =  1,M 

FILAIR(I)  =  FILNAM(I) 
FILDAT(I)  =  FILNAM(I) 
FILMIN(I)  =  FILNAM(I) 
FILPAT(I)  =  FILNAM(I) 

FILL  OUT  THE  REST  OF  THE  FILE  NAMES 

AIRWAY  FILE 


FILAIR(M+1) 

= 

ICHA 

FILAIR(M+2) 

= 

ICHI 

FILAIR(M+3) 

= 

ICHR 

FILAIR(M+4) 

= 

I  DOT 

FILAIR(M+5) 

= 

ICHA 

FILAIR(M+6) 

s 

ICHS 

FILAIR(M+7) 

= 

ICHC 

FILAIR(M+8) 

= 

0 

BINARY  FILE 

FILDAT(M+1) 

a 

ICHD 

FILDAT(M+2) 

= 

ICHA 
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0172 
0173 
0174 
0175 
0176 
0177 
0178 
0179 
0180 
0181 
0182 
0183 
0184 
0185 
0186 
0187 
0188 
0189 
0190 
0191 
0192 
0193 
0194 
0195 
0196 
0197 
0198 
0199 
0200 
0201 
0202 
0203 
0204 
0205 
0206 
0207 
0208 
0209 
0210 
0211 
0212 
0213 
0214 
0215 
0216 
0217 
0218 
0219 
0220 
0221 
0222 
0223 
0224 
0225 
0226 
0227 
0228 
0229 


C 
C 
C 

137 

C 
C 
C 


3100 
3001 

1001 


FILDAT(M+3) 

= 

ICHT 

FILDAT(M+4) 

= 

I  DOT 

FILDAT(M+5) 

= 

ICHD 

FILDAT(M+6) 

= 

ICHA 

FILDAT(M+7) 

= 

ICHT 

FILDAT(M+8) 

= 

0 

MINER  FILE 


FILMIN( 
FILMIN( 
FILMIN( 
FILMIN( 
FILMIN( 
FILMIN( 
FILMIN( 
FILMIN( 


M+l) 
M+2) 
M+3) 
M+4) 
M+5) 
M+6) 
M+7) 
M+8) 


ICHM 
ICHI 
ICHN 
I  DOT 
ICHA 
ICHS 
ICHC 
0 


ESCAPE  PATH  FILE 


FILPAT(M+1)  = 

ICHP 

FILPAT(M+2)  = 

ICHA 

FILPAT(M+3)  = 

ICHT 

FILPAT(M+4)  = 

IDOT 

FILPAT(M+5)  = 

ICHA 

FILPAT(M+6)  = 

ICHS 

FILPAT(M+7)  = 

ICHC 

FILPAT(M+8)  = 

0 

OPEN  FILE  FOR 

BINARY  OUTPUT 

DATA 


OPEN  (   UNIT=FDAT, 

1        NAME  =  FILDAT, 

1        RECORDSIZE  =  8, 

1        TYPE  =  'UNKNOWN' , 

1        FORM  =  'UNFORMATTED', 

1        ACCESS  =  'DIRECT' ) 

NULL  BINARY  FILE 

DO  137  1=1,300 

WRITE  (FDAT' I ) IFILL4 , IFILL4 , IFILL4 , IFILL4 , 

1        IFILL4,IFILL4,IFILL4,IFILL4 

OPEN  FILE  FOR  AIRWAY  ASCII  INPUT  DATA 

OPEN(UNIT  =  FAIR, 

1  NAME  =  FILAIR, 

2  TYPE  =  'OLD' ) 
WRITE  (5,3100) 
FORMAT ( '  ' ) 
WRITE  (5,3001) 

FORMAT ( *   PROCESSING  AIRWAY  FILE') 

DO  1001  I  =  1,100 

READ (FAIR,*) IAIR, I JUN1 , I JUN2 , HEIGHT , LENGTH 

IF  (IAIR.EQ.999)GO  TO  200 

WRITE( FDAT' 100+IAIR) 

1      IAIRrIJUNl,FILL4,IJUN2, FILL4,IFILL4, HEIGHT, LENGTH, I FILL2 
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0230  C 

0231  C        CLOSE  AIRWAY  FILE 

0232  C 

0233  200      CLOSE  (UNIT  =  FAIR) 

0234  C 

0235  C        OPEN  FILE  FOR  MINER  ASCII  INPUT  DATA 

0236  C 

0237  OPEN(UNIT  =  FMIN, 

0238  1  NAME  =  FILMIN, 

0239  2  TYPE  =  'OLD* ) 

0240  2000     CONTINUE 

0241  WRITE  (5,3100) 

0242  WRITE  (5,3002) 

0243  3002     FORMAT('   PROCESSING  MINER  FILE') 

0244  DO  2001  I  =  1,100 

0  245  READ(FMIN,*)IMIN,INILOC,IESPAT,VD02M 

0246  IF  (IMIN.EQ.999)GO  TO  202 

0247  WRITE(FDAT' IMIN) 

0248  1        IMIN, FILL, INILOC, FILL, IESPAT,VD02M , IFILL4 , FILL, IFILL2 

0249  READ( FMIN , * ) IMIN ,WMASS , IRJUNl , RT1 , IRJUN2 , RT2 , IRJUN3 , RT3 

0250  IF  (IMIN.EQ.999)GO  TO  202 

0251  2001     WRITE(FDAT'50  +  IMIN) 

0252  1        IMIN, WMASS, IRJUNl, RT1 , IRJUN2 ,RT2 , IRJUN3 , RT3 , IFILL2 

0253  C 

0254  C        CLOSE  MINER  FILE 

0255  C 

0256  202      CLOSE  (UNIT  =  FMIN) 

0257  C 

0  258  C        OPEN  FILE  FOR  PATH  ASCII  INPUT  DATA 

0259  C 

0260  OPEN(UNIT  =  FPAT, 

0261  1  NAME  =  FILPAT, 

0262  2  TYPE  =  'OLD' ) 

0263  WRITE  (5,3100) 

0264  WRITE  (5,3003) 

0265  3003     FORMAT('   PROCESSING  PATH  FILE') 

0266  DO  4001  I  =  1,100 

0267  READ( FPAT, * ) IPAT, IAIR1 , IAIR2 , IAIR3 , IAIR4 , IAIR5 , 

0268  1        IAIR6,IAIR7,NXTPAT 

0269  IF  (IPAT.EQ.999)GO  TO  203 

0270  4001     WRITE (FDAT' 200  +  IPAT) 

0271  1        IPAT,IAIR1,IAIR2,IAIR3,IAIR4,IAIR5,IAIR6,IAIR7,NXTPAT 

0272  C 

0273  C        CLOSE  FILES 

0274  C 

0275  203      CLOSE  (UNIT  =  FPAT) 

0276  CLOSE  (UNIT  =  FDAT) 

0277  WRITE  (5,3100) 

0278  STOP  '  MINE  DATA  FILE  TRANSLATED  ' 

0279  END 
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APPENDIX  C— LISTING  OF  RESCUE. FOR 

!       THIS  DIALOGUE  GENERATED  THE  OUTPUT  SHOWN  IN 
!       APPENDIX  E 
$  RUN  RESCUE 

OUTPUT  DEVICE 

CURRENT  VALUE  OF  LP  IS      5 

NEW  LP?   RETURN  FOR  NO  CHANGE  3 

CURRENT  VALUE  OF  LP  IS      3 

ARE  REST  JUNCTIONS  TO  BE  INCLUDED? 
JFLG=1   :   OMIT  REST  JUNCTIONS 
JFLG=0   :   INCLUDE  REST  JUNCTIONS 

CURRENT  VALUE  OF  JFLG  IS      0 

NEW  JFLG?   RETURN  FOR  NO  CHANGE  1 

CURRENT  VALUE  OF  JFLG  IS      1 

ARE  MINERS  ALLOWED  TO  RUN  WHEN  HEIGHT  >  84"? 
IWORR  =  0  :  MINERS  ARE  TO  WALK 
IWORR  =  1  :  MINERS  ARE  TO  RUN 

CURRENT  VALUE  OF  IWORR  IS      0 

NEW  IWORR?   RETURN  FOR  NO  CHANGE  0 

CURRENT  VALUE  OF  IWORR  IS      0 

ENTER  PREFIX  FOR  DATA  FILE  NAMES:   MN1 

VELOCITY  CONSTANTS 

THE  RELATIONSHIP  BETWEEN  VELOCITY  AND  SEAM  HEIGHT  IS 

GIVEN  BY: 

V=A+B*H+C*H**2+D*H**3 

WHERE  H  IS  IN  INCHES,  AND  V  IS  IN  FEET/MIN 

A  =  252.18 

NEW  A?  RETURN  FOR  NO  CHANGE 
A  =   252.2 
B  =  -16.139 

NEW  B?  RETURN  FOR  NO  CHANGE 
B  =  -16.14 
C  =  0.408 
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NEW  C?  RETURN  FOR  NO  CHANGE 
C  «  0.4080 
D  »  -0.00244 

NEW  D?  RETURN  FOR  NO  CHANGE 
D  *  -.2440E-02 

0-2  CONSUMPTION  EQUATION  FOR  CRAWLING 

VDOTO-2  =  0.00  *  (  M  **  0.00  )  *  VP  **  0 
+  48.60  *  (  M  **  0.00  )  *  VP  **  1 
+   0.00  *  (  M  **   0.00  )  *  VP  **  % 

WHERE  VDOTO-2  IS  IN  ML(STP)/KG/MIN, 

M  IS  IN  KG,  AND  VP  IS  IN  M/SEC 

DO  YOU  WANT  TO  CHANGE  THIS  EXPRESSION  ?  (Y  OR  N) 

0-2  CONSUMPTION  EQUATION  FOR  DUCKWALK 

VDOTO-2  =   0.00  *  (  M  **   0.00  ) 

+  39.60  *  (  M  **   0.00  ) 

+   0.00  *  (  M  **   0.00  ) 

WHERE  VDOTO-2  IS  IN  ML( STP)/KG/MIN, 

M  IS  IN  KG,  AND  VP  IS  IN  M/SEC 

DO  YOU  WANT  TO  CHANGE  THIS  EXPRESSION  ?  (Y  OR  N) 

0-2  CONSUMPTION  EQUATION  FOR  WALKING 

VDOTO-2  =  7.06  *  (  M  **  0.00  )  *  VP  **  0 
+  0.00  *  (  M  **  0.00  )  *  VP  **  1 
+   3.60  *  (  M  **   0.00  )  *  VP  **  2 

WHERE  VDOTO-2  IS  IN  ML(STP)/KG/MIN, 

M  IS  IN  KG,  AND  VP  IS  IN  M/SEC 

DO  YOU  WANT  TO  CHANGE  THIS  EXPRESSION  ?  (Y  OR  N) 

0-2  CONSUMPTION  EQUATION  FOR  RUNNING 

VDOTO-2  =  4.50  *  (  M  **  0.00  )  *  VP  **  0 
+  12.00  *  (  M  **  0.00  )  *  VP  **  1 
+   0.00  *  (  M  **   0.00  )  *  VP  **  2 

WHERE  VDOTO-2  IS  IN  ML( STP) /KG/MIN , 

M  IS  IN  KG,  AND  VP  IS  IN  M/SEC 

DO  YOU  WANT  TO  CHANGE  THIS  EXPRESSION  ?  (Y  OR  N) 


* 

VP 

**    0 

* 

VP 

**    ^ 

* 

VP 

**    2 
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"S"  TO  STOP,  RETURN  TO  CONTINUE 

"S"  TO  STOP,  RETURN  TO  CONTINUE 

"S"  TO  STOP,  RETURN  TO  CONTINUE 

"S"  TO  STOP,  RETURN  TO  CONTINUE 

"S"  TO  STOP,  RETURN  TO  CONTINUE 
FORTRAN  STOP 
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0001 

C 

0002 

C 

0003 

C 

0004 

C 

0005 

C 

0006 

c 

0007 

c 

0008 

c 

0009 

c 

0010 

c 

0011 

c 

0012 

c 

0013 

c 

0014 

c 

0015 

c 

0016 

c 

0017 

c 

0018 

c 

0019 

c 

0020 

c 

0021 

c 

0022 

c 

0023 

c 

0024 

c 

0025 

c 

0026 

c 

0027 

c 

0028 

c 

0029 

c 

0030 

c 

0031 

c 

0032 

c 

0033 

c 

0034 

c 

0035 

c 

0036 

c 

0037 

c 

0038 

c 

0039 

c 

0040 

c 

0041 

c 

0042 

c 

0043 

c 

0044 

c 

0045 

c 

0046 

c 

0047 

c 

0048 

c 

0049 

c 

0050 

c 

0051 

c 

0052 

c 

0053 

c 

0054 

c 

0055 

c 

0056 

c 

0057 

c 

THIS  PROGRAM  CALCULATES  WHETHER  OR  NOT  A  MINER  WOULD 
BE  ABLE  TO  MAKE  IT  TO  SAFTEY  FOR  A  GIVEN  ESCAPE  PATH, 
OXYGEN  USAGE,  AND  EMERGENCY  SUPPLY  OF  OXYGEN. 

SET  UP  FOR: 

50  MINERS 

ITEM  VARIABLE 

FIRST  50  RECORDS  OF  %%%DAT.DAT 


STARTING  JUNCTION 
ESCAPE  PATH 
0-2  RATE 

SECOND  50  RECORDS  OF  %%%DAT.DAT 

MASS 

RESTING  JUNCTION 
RESTING  TIME 
RESTING  JUNCTION 
RESTING  TIME 
RESTING  JUNCTION 
RESTING  TIME 

100  AIRWAYS 

ITEM 

JUNCTION  1 

JUNCTION  2 

JUNCTION  3 

JUNCTION  4 

JUNCTION  5 

JUNCTION  6 
HEIGHT 

50  EXCAPE  PATHS 

ITEM 

STARTING  JUNCTION 
ENDING  JUNCTION 
NEXT  PATH 


JUNSTR 
ESCPAT 
VDOT 


MASS 

RESJUN(I,1) 

RESTIM(I,1) 

RESJUN(If2) 

RESTIM(I,2) 

RESJUN(I,3) 

RESTIM(I,3) 


VARIABLE 

AIRJUN(I,1) 
AIRJUN(I,2) 
AIRJUN(I,3) 
AIRJUN(I,4) 
AIRJUN(I,5) 
AIRJUN(I,6) 
AIRHGT 


VARIABLE 

ESCSTR 
ESCEND 
NXTPAT 


THIS  DATA  IS  STORED  IN  FILE  %%%DAT.DAT  (THE  FIRST  3  CHARACTERS  OF 
THE  NAME  INDICATE  TO  WHICH  MINE  THE  FILE  REFERS.) 

THEY  ARE  TRANSLATED  INTO  THIS  BINARY  FILE  BY  THE  PROGRAM  TRANS. 
THE  EDITOR  IS  USED  TO  CREATE  THE  ASCII  VERSION  OF  THE  DATA  FILES. 
THESE  FILES  ARE  THEN  RUN  THROUGH  TRANS  AND  A  BINARY  VERSION  IS 
WRITTEN.   THERE  ARE  3  ASCII  FILES  USED  BY  THIS  PROGRAM. 

THESE  FILES  ARE  (THE  FIRST  3  CHARACTERS  OF  THE  NAME  INDICATE 
TO  WHICH  MINE  THE  FILE  REFERS,  MINE  1  IN  THIS  EXAMPLE.) 

MN1AIR.ASC 


0058 

C 

0059 

C 

0060 

C 

0061 

C 

0062 

C 

0063 

C 

0064 

C 

0065 

C 

0066 

c 

0067 

c 

0068 

c 

0069 

c 

0070 

c 

0071 

c 

0072 

c 

0073 

c 

0074 

c 

0075 

c 

0076 

c 

0077 

c 

0078 

c 

0079 

c 

0080 

c 

0081 

c 

0082 

0083 

c 

0084 

c 

0085 

c 

0086 

0087 

0088 

0089 

0090 

0091 

0092 

c 

0093 

c 

0094 

c 

0095 

0096 

0097 

0098 

c 

0099 

c 

0100 

c 

0101 

0102 

0103 

c 

0104 

c 

0105 

c 

0106 

c 

0107 

c 

0108 

c 

0109 

0110 

0111 

c 

0112 

c 

0113 

c 

0114 
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THIS  IS  AN  AIRWAY  LIST  FOR  MINE  1. 

MN1MIN.ASC 
THIS  IS  A  MINER  LIST  FOR  MINE  1. 

MN1PAT.ASC 

THIS  IS  A  ESCAPE  PATH  LIST  FOR  MINE  1. 

THESE  ASCII  FILES  ARE  TRANSLATED  INTO  A  DATA  BASE  FILE  MN1DAT.DAT. 
THE  .DAT  FILE  HAS  HE  FOLLOWING  STRUCTURE: 

EACH  RECORD  IS  3  2  BYTES  (8  LONG  WORDS)  LONG. 

50  RECORDS  FOR  LINE  1  DATA  FROM  MN1MIN.ASC 

50  RECORDS  FOR  LINE  2  DATA  FROM  MN1MIN.ASC 

100  RECORDS  FOR  DATA  FROM  MN1AIR.ASC 

100  RECORDS  FOR  DATA  FROM  MNlPAT.ASC 

PROGRAM  RESCUE 

MINER  RESERVES 

REAL     *  4  MASS(0:50) 

REAL     *  4  VDOT(0:50) 

INTEGER  *  4  JUNSTR(0:50) 

INTEGER  *  4  ESCPAT(0:50) 

INTEGER  *  4  RES JUN( 0 : 50 , 3 ) 

REAL     *  4  RESTIM(0:50,3) 

AIRWAY  RESERVES 

INTEGER  *  4  AIRJUN ( 0 : 100 , 6 ) 
REAL     *  4  AIRHGT( 0:100) 
REAL     *  4  AIRLEN(0 :100) 

PATH  RESERVES 

INTEGER  *  4  ESCAIR( 0 : 50 , 7 ) 
INTEGER  *  4  NXTPAT(0:50) 

INTERIUM  CALCULATIONS 

ALLOW  21  AIRWAYS  AND  3  REST  JUNCTIONS 
TO  REACH  FRESH  AIR 

REAL  *4  TIME(24) 
REAL  *4  U(24) 

FILE  NAME 

BYTE  FILNAM(12) 
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0115 

0116 

C 

0117 

C 

0118 

C 

0119 

0120 

0121 

0122 

C 

0123 

C 

0124 

C 

0125 

0126 

C 

0127 

C 

0128 

C 

0129 

0130 

C 

0131 

C 

0132 

C 

0133 

0134 

C 

0135 

C 

0136 

c 

0137 

0138 

c 

0139 

c 

0140 

c 

0141 

0142 

c 

0143 

0144 

c 

0145 

0146 

c 

0147 

c 

0148 

c 

0149 

0150 

0151 

0152 

0153 

c 

0154 

0155 

c 

0156 

0157 

0158 

0159 

0160 

0161 

c 

0162 

0163 

0164 

c 

0165 

c 

0166 

c 

0167 

0168 

0169 

0170 

0171 

c 

BYTE  FILDATQ2) 

VARIABLES  TO  INDEX  AIRWAYS,  MINERS,  AND  PATHS 

INTEGER  *2  IAIR 
INTEGER  *2  IMIN 
INTEGER  *2  IPAT 

VARIABLES  FOR  LUNS 

INTEGER  *2  FDAT 

TWO  BYTE  FILL  CHARACTER 

INTEGER  *2  IFILL2 

FOUR  BYTE  INTEGER  FILL  CHARACTER 

INTEGER  *4  IFILL4 

FOUR  BYTE  REAL  FILL  CHARACTER 

REAL  *4  FILL 

READ  BUFFER  RESERVES 

REAL  *4  RELBUF(4) 

INTEGER  *4  I4BUF(7) 

INTEGER  *2  ISBUF,IEBUF 

DATA  REDUCTION  CONSTANT  RESERVES 


REAL  *4  A 
REAL  *4  B 
REAL  *4  C 
REAL  *4  D 

BYTE  IANS(5) 


BYTE 
BYTE 
BYTE 
BYTE 
BYTE 

REAL 
REAL 


ITEXT(4,8) 

ITEXTK8) 

ITEXT2(8) 

ITEXT3(8) 

ITEXT4(8) 

*4  ALPHA(3,4) 
*4  LAMDA(3,4) 


! CONVERSION 
1FROM  HEIGHT 
!TO 
[VELOCITY 

! TERMINAL  RESPONSE 

•TEXT  TO  OUTPUT 
!  INPUT  TEXT 
! INPUT  TEXT 
!  INPUT  TEXT 
■INPUT  TEXT 

1COEFFICIENT  MATRIX  FOR  0-2  USE 
! EXPONENT  MATRIX  FOR  0-2  USE 


SET  UP  TEXTS  FOR  LATER  PRINTOUT 

DATA  ITEXT1  / ' C ' , ' R ' , ' A' , ' W , ' L ' , ' I ' , * N ' , ' G ' / 

DATA  ITEXT2  / ' D' , ' U ' , ' C ' , ' K ' , ' W' , ' A' , ' L * , * K' / 

DATA  ITEXT3  /'W  , ' A*  , ' L '  , ' K '  , ' I  *  , ' N '  , ' G *  , '  '/ 

DATA  ITEXT4  /* R' , ' U ' , ' N ' , ' N ' , ' I ' , ' N ' , ' G ' , '  '/ 
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0172 

C 

GET  TEXT  INTO  ITEXT  IN  A  WAY  THAT  MAKES  IT  EA 

0173 

C 

AND  STILL  RETAINS  A  READABLE  INPUT  FORMAT 

0174 

C 

0175 

DO  I  =  1,8 

0176 

ITEXT (1,1)  =  ITEXTl(I) 

0177 

ITEXT(2,I)  =  ITEXT2(I) 

0178 

ITEXT (3, I)  =  ITEXT3(I) 

0179 

ITEXT (4, I)  =  ITEXT4(I) 

0180 

END  DO 

0181 

C 

0182 

C 

NULL  FILL  CHARACTERS 

0183 

C 

0184 

LP=5 

0185 

IFILL2  =  0 

0186 

IFILL4  =  0 

0187 

FILL  =  0. 

0188 

C 

0189 

C 

5  MPH  IN  M/SEC 

0190 

C 

0191 

C 

MPH       88  FT/SEC      12  IN      0.0254  M 

0192 
0193 

C 
C 

* * * 

60  MPH          FT          IN 

0194 

C 

0195 

V5MPH  =  5.  *  88.  *  12.  *  0.0254  /60. 

0196 

C 

0197 

C 

FIND  OUT  ABOUT  OUTPUT  DEVICE 

0198 

C 

0199 

WRITE(5,2000) 

0200 

WRITE (5, 5) 

0201 

5 

FORMAT ('  OUTPUT  DEVICE') 

0202 

WRITE (5, 7)   LP 

0203 

7 

FORMAT('  CURRENT  VALUE  OF  LP  IS  ' ,1X,I4) 

0204 

CALL  PRMPTF('  NEW  LP?   RETURN  FOR  NO  CHANGE  ' 

0205 

IF(NF.NE.O)  LP=XPNEW 

0206 

WRITE(5,2000) 

0207 

WRITE (5, 7)  LP 

0208 

C 

0209 

C 

FIND  OUT  IF  WE  ARE  TO  INCLUDE  REST  JUNCTIONS 

0210 

C 

JFLG=1  :  OMIT  REST  JUNCTIONS 

0211 

C 

JFLG=0  :  INCLUDE  REST  JUNCTIONS 

0212 

c 

0213 

WRITE(5,2000) 

0214 

JFLG=0 

0215 

WRITE(5,19) 

0216 

19 

FORMAT ('  ARE  REST  JUNCTIONS  TO  BE  INCLUDED?  ' 

0217 

1'       JFLG=1   :   OMIT  REST  JUNCTIONS'/ 

0218 

2'       JFLG=0   :   INCLUDE  REST  JUNCTIONS') 

0219 

WRITE(5,2000) 

0220 

WRITE(5,22)  JFLG 

0221 

WRITE(5,20T)0) 

0222 

22 

FORMAT ( '  CURRENT  VALUE  OF  JFLG  IS  ',1X,I4) 

0223 

CALL  PRMPTF('  NEW  JFLG?   RETURN  FOR  NO  CHANGE 

0224 

IF(NF.NE.O)  JFLG=XJFLG 

0225 

WRITE(5,2000) 

0226 

WRITE(5,22)  JFLG 

0  227 

WRITE(LP,22)  JFLG 

0228 

WRITE(LP,2000) 

,XPNEW,NF) 


,XJFLG,NF) 
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0229 

C 

0230 

C 

0231 

C 

0232 

C 

0233 

C 

0234 

0235 

0236 

0237 

23 

0238 

0239 

0240 

0241 

0242 

0243 

24 

0244 

0245 

0246 

0247 

0248 

0249 

0250 

0251 

C 

0252 

C 

0253 

C 

0254 

0255 

200 

0256 

C 

0257 

c 

0258 

c 

0259 

0260 

c 

0261 

c 

0262 

c 

0263 

1 

0264 

c 

0265 

c 

0266 

c 

0267 

0268 

100 

0269 

c 

0270 

c 

0271 

c 

0272 

c 

0273 

c 

0274 

0275 

0276 

0277 

0278 

0279 

0280 

0281 

0282 

c 

0283 

c 

0284 

c 

0285 

FIND  OUT  IF  MINERS  ARE  ALLOWED  TO  RUN  WHEN  HEIGHT  >  84" 
IWORR  =  0  :  MINERS  ARE  TO  WALK 
IWORR  =  1  :  MINERS  ARE  TO  RUN 

WRITE(5,2000) 

IWORR  =  0 

WRITE(5,23) 

FORMAT ( '  ARE  MINERS  ALLOWED  TO  RUN  WHEN  HEIGHT  >  84"?  '/ 

1'       IWORR  =  0  :  MINERS  ARE  TO  WALK  '/ 

2'       IWORR  =  1  :  MINERS  ARE  TO  RUN   ') 

WRITE(5,2000) 

WRITE (5, 24)  IWORR 

WRITE (5, 2000) 

FORMAT('  CURRENT  VALUE  OF  IWORR  IS  ' ,1X,I4) 

CALL  PRMPTF('  NEW  IWORR?   RETURN  FOR  NO  CHANGE  ',XIWORR,NF) 

IF(NF.NE.O)  IWORR  =  XIWORR 

WRITE(5,2000) 

WRITE (5, 24)  IWORR 

WRITE(LP,24)  IWORR 

WRITE(LP,2000) 

WRITE  (5,2000) 

NULL  CERTAIN  VARIABLES 

DO  200  I  =  1,50 
JUNSTR(I)  =  0 

SET  UP  LUNS 

FDAT  =  8  !LUN  FOR  BINARY  FILE 

GET  FILE  NAME  FOR  DATA  INPUT 

CALL  PRMPTA( '  ENTER  PREFIX  FOR  DATA  FILE  NAMES:   ',FILNAM,M) 

TRANSFER  PREFIX  OF  FILE  NAMES  INTO  REQUIRED  VARIABLES 

DO  100  I  =  1,M 
FILDAT(I)  =  FILNAM(I) 

FILL  OUT  THE  REST  OF  THE  FILE  NAMES 

BINARY  FILE 

FILDAT(M+1)  =  'D' 

FILDAT(M+2)  =  'A' 

FILDAT(M+3)  =  'T' 

FILDAT(M+4)  =  ' . * 

FILDAT(M+5)  =  'D' 

FILDAT(M+6)  =  'A' 

FILDAT(M+7)  =  »T' 

FILDAT(M+8)  =  0 

OPEN  FILE  FOR  BINARY  INPUT  DATA 

OPEN  (   UNIT=FDAT, 


23 


0286 

1        NAME  =  FILDAT, 

0287 

1        RECOPDSIZE  =  8, 

0288 

1        TYPE  =  'OLD' , 

0289 

1        FORM  =  'UNFORMATTED', 

0290 

1        ACCESS  =  'DIRECT' ) 

0291 

C 

0292 

C 

READ  DATA  FROM  BINARY  FILE 

0293 

C 

0294 

C 

AIRWAY  DATA 

0295 

C 

0296 

C 

JUNCTION  1 

0297 

C 

JUNCTION  2 

0298 

C 

HEIGHT 

0299 

c 

LENGTH 

0300 

c 

0301 

DO  10  I  =  1,100 

0302 

c 

0303 

c 

READ  DATA 

0304 

c 

0305 

READ(FDAT'100  +  DISBUF, 

0306 

1        I4BUF(1) ,RELBUF(1) , 

0307 

2        I4BUF(2) ,RELBUF(2) , 

0308 

3        I4BUF(3) ,RELBUF(3) , 

0309 

5        RELBUF(4) ,IEBUF 

0310 

c 

0311 

c 

TRANSFER  TO  VARIABLES 

0312 

c 

0313 

AIRHGT  (ISBUF)  =  RELBUF(3) 

0314 

AIRLEN  (ISBUF)  =  RELBUF(4) 

0315 

DO  10  11  =  1,2 

0316 

AIRJUN(I,I1)  =  I4BUF(I1) 

0317 

10 

CONTINUE 

0318 

c 

0319 

c 

READ  MINER  DATA  (FIRST  PART) 

0320 

c 

0321 

c 

FIRST  50  RECORDS  OF  %%%DAT.D£ 

0322 

c 

0323 

c 

STARTING  JUNCTION 

0324 

c 

ESCAPE  PATH 

0325 

c 

0-2  RATE 

0326 

c 

0327 

DO  30  I  =  1,50 

0328 

c 

0329 

c 

READ  DATA 

0330 

c 

0331 

READ(FDAT'  DISBUF, 

0332 

1        RELBUF(l) ,I4BUF(1) , 

0333 

2        RELBUF(2) ,I4BUF(2) , 

0334 

3        RELBUF(3)  ,I4BUF(3)  , 

0335 

4        RELBUF(4), 

0336 

5        IEBUF 

0337 

c 

0338 

c 

TRANSFER  TO  VARIABLES 

0339 

c 

0340 

JUNSTR  (ISBUF)  =  I4BUF(1) 

0341 

ESCPAT  (ISBUF)  =  I4BUF(2) 

0342 

VDOT  (ISBUF)  =  RELBUF(3) 

AIRJUN(I,1) 
AIRJUN(I,2) 
AIRHGT 
AIRLEN 


JUNSTR 
ESCPAT 
VDOT 
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0343 

30 

0344 

C 

0345 

C 

0346 

C 

0347 

C 

0348 

C 

0349 

C 

0350 

C 

0351 

C 

0352 

C 

0353 

C 

0354 

C 

0355 

c 

0356 

c 

0357 

0358 

c 

0359 

c 

0360 

c 

0361 

0362 

0363 

0364 

0365 

0366 

0367 

c 

0368 

c 

0369 

c 

0370 

0371 

0372 

0373 

0374 

0375 

0376 

0377 

40 

0378 

C 

0379 

c 

0380 

c 

0381 

c 

0382 

c 

0383 

c 

0384 

c 

0385 

c 

0386 

c 

0387 

c 

0388 

c 

0389 

c 

0390 

0391 

c 

0392 

c 

0393 

c 

0394 

0395 

0396 

0397 

0398 

0399 

CONTINUE 

READ  MINER  DATA  (SECOND  PART) 

SECOND  50  RECORDS  OF  %%%DAT.DAT 

MASS 

RESTING  JUNCTION 
RESTING  TIME 
RESTING  JUNCTION 
RESTING  TIME 
RESTING  JUNCTION 
RESTING  TIME 

DO  40  I  =  1,50 

READ  DATA 

READ(FDAT'50  +  I) ISBUF, 

1  RELBUF(l) ,I4BUF(1) , 

2  RELBUF(2) ,I4BUF(2) , 

3  RELBUF(3) ,I4BUF(3) , 

4  RELBUF(4), 

5  IEBUF 

TRANSFER  TO  VARIABLES 

MASS  (ISBUF)  =  RELBUF(l) 

RESJUN  ( ISBUF,  1)  =  I4BUFQ) 

RESTIM  ( ISBUF, 1)  =  RELBUF(2) 

RESJUN  (ISBUF, 2)  =  I4BUF(2) 

RESTIM  (ISBUF, 2)  =  RELBUF(3) 

RESJUN  (ISBUF, 3)  =  I4BUF(3) 

RESTIM  (ISBUF, 3)  =  RELBUF(4) 
CONTINUE 

READ  ESCAPE  PATH  DATA 


DO  50  I  =  1,50 

READ  DATA 

READ(FDAT'200  +  I) ISBUF, 


MASS 

RESJUN (1,1) 
RESTIM(I,1) 
RESJUN (I, 2) 
RESTIM(I,2) 
RESJUN (I, 3) 
RESTIM(I,3) 


ESCAPE 

AIRWAY 

1 

ESCAIRd 

,1) 

ESCAPE 

AIRWAY 

2 

ESCAIRd, 

,2) 

ESCAPE 

AIRWAY 

3 

ESCAIRd 

r3) 

ESCAPE 

AIRWAY 

4 

ESCAIRd, 

,4) 

ESCAPE 

AIRWAY 

5 

ESCAIRd 

r5) 

ESCAPE 

AIRWAY 

6 

ESCAIRd, 

r6) 

ESCAPE 

AIRWAY 

7 

ESCAIRd 

r7) 

NEXT  PATH 

NXTPAT 

I4BUF(1) ,I4BUF(2) , 
I4BUF(3) ,I4BUF(4) , 
I4BUF(5) ,I4BUF(6) , 
I4BUF(7) , 
IEBUF 
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0400  C 

0401  C        TRANSFER  TO  VARIABLES 

0402  C 

0403  DO  50  II  =  1,7 

0404  ESCAIR(ISBUF,I1)  =  I4BUF(I1) 

0405  NXTPAT  (ISBUF)  =  IEBUF 

0406  50       CONTINUE 

0407  C 

0408  C        VELOCITY  CONSTANTS 

0409  C        THE  RELATIONSHIP  BETWEEN  VELOCITY  AND  SEAM  HEIGHT  IS 

0410  C        GIVEN  BY: 

0411  C  V=A+B*H+C*H**2+D*H**3 

0412  C 

0413  C        WHERE  H  IS  IN  INCHES,  AND  V  IS  IN  FEET/MIN 

0414  C 

0415  WRITE  (5,2000) 

0416  WRITE  (5,2007) 

0417  2007     FORMAT('  VELOCITY  CONSTANTS ' ) 

0418  WRITE  (5,2008) 

0419  2008     FORMAT ( '  THE  RELATIONSHIP  BETWEEN  VELOCITY  AND  SEAM  HEIGHT  IS  ') 

0420  WRITE  (5,2009) 

0421  2009     FORMAT*'  GIVEN  BY:') 

0422  WRITE  (5,2010) 

0423  2010     FORMAT ( '         V=A+B*H+C*H**2+D*H**3') 

0424  WRITE  (5,2000) 

0425  WRITE  (5,2011) 

0426  20.11     FORMAT('  WHERE  H  IS  IN  INCHES,  AND  V  IS  IN  FEET/MIN') 

0427  WRITE  (5,2000) 

0428  C 

0429  A  =  252.18 

0430  B  =  -16.139 

0431  C  =  0.408 

0432  D  =  -0.00244 

0433  C 

0434  WRITE  (5,2012) 

0435  2012     FORMAT ( '  A  =  252.18') 

0436  CALL  PRMPTF('  NEW  A?  RETURN  FOR  NO  CHANGE  ',ANEW,NF) 

0437  IF  (NF.NE.O)A  =  ANEW 

0438  WRITE  (LP, 2013)  A 

0439  WRITE  (5,2000) 

0440  WRITE  (5,2013)A 

0441  2013     FORMAT ( '  A  =  ',G10.4) 

0442  WRITE  (5,2000) 

0443  WRITE  (5,2014) 

0444  2014     FORMATC  B  =  -16.139') 

0445  CALL  PRMPTF('  NEW  B?  RETURN  FOR  NO  CHANGE  ',BNEW,NF) 

0446  IF  (NF.NE.O)B  =  BNEW 

0447  WRITE  (5,2000) 

0448  WRITE  (LP,2015)B 

0449  WRITE  (5,2015)B 

0450  2015     FORMAT('  B  =  *,G10.4) 

0451  WRITE  (5,2000) 

0452  WRITE  (5,2016) 

0453  2016     FORMAT ( '  C  =  0.408') 

0454  CALL  PRMPTF('  NEW  C?  RETURN  FOR  NO  CHANGE  ',CNEW,NF) 

0455  IF  (NF.NE.O)C  =  CNEW 

0456  WRITE  (5,2000) 
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0457 
0458 
0459 
0460 
0461 
0462 
0463 
0464 
0465 
0466 
0467 
0468 
0469 
0470 
0471 
0472 
0473 
0474 
0475 
0476 
0477 
0478 
0479 
0480 
0481 
0482 
0483 
0484 
0485 
0486 
0487 
0488 
0489 
0490 
0491 
0492 
0493 
0494 
0495 
0496 
0497 
0498 
0499 
0500 
0501 
0502 
0503 
0504 
0505 
0506 
0507 
0508 
0509 
0510 
0511 
0512 
0513 


2017 


2018 


2019 

C 
C 
C 
C 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


WRITE  (LP, 2017)  C 

WRITE  (5,2017)C 

FORMAT( '  C  =  ' ,G10.4) 

WRITE  (5,2000) 

WRITE  (5,2018) 

FORMAT(  *  D  =  -0.00244'  ) 

CALL  PRMPTF('  NEW  D?  RETURN  FOR  NO  CHANGE  ',DNEW,NF) 

IF  (NF.NE.O)D  =  DNEW 

WRITE  (5,2000) 

WRITE(LP,2019)  D 

WRITE  (5,2019)D 

FORMAT ( '  D  =  ' ,G10.4) 

WRITE  (5,2000) 

0-2  CONSUMPTION  CONSTANTS 

THE  RELATIONSHIP  BETWEEN  0-2  CONSUMPTION,  MASS,  AND  VELOCITY 

IS  GIVEN  BY: 

VDOT02  =  ALPHA (1, J) 

+  ALPHA(2,J) 

+  ALPHA(3,J) 


*  (  M  **  LAMDA(1,J) )  * 

*  (  M  **  LAMDA(2,J) )  * 

*  (  M  **  LAMDA(3,J) )  * 


VP 
VP 
VP 


** 
** 
** 


WHERE  M  IS  IN  KG,  AND  VDOT02  IS  IN  ML  O  -  2/KG/MIN 
VP  IS  IN  M/SEC,  AND  J  INDICATES  AS  FOLLOWS: 

1  =>  CRAWL 

2  =>  DUCKWALK 

3  =>  WALK 

4  =>  RUN 


CRAWL 


ALPHA(1,1) 

= 

0. 

LAMDA(1,1) 

= 

0. 

ALPHA(2,1) 

= 

48.6 

LAMDA(2,1) 

= 

0. 

ALPHA(3,1) 

= 

0. 

LAMDA(3,1) 

= 

0. 

DUCKWALK 

ALPHA(1,2) 

= 

0. 

LAMDA (1,2) 

= 

0. 

ALPHA(2,2) 

= 

39.6 

LAMDA(2,2) 

= 

0. 

ALPHA(3,2) 

= 

0. 

LAMDA(3,2) 

= 

0. 

WALKING 

ALPHA (1,3) 

= 

7.06 

LAMDA (1,3) 

= 

0. 

ALPHA(2,3) 

= 

0. 

LAMDA (2, 3) 

= 

0. 

ALPHA(3,3) 

= 

3.6 

LAMDA (3, 3) 

= 

0. 

RUNNING 

!V 

** 

0 

!V 

** 

0 

IV 

** 

1 

!V 

*  * 

1 

•V 

*  * 

2 

!V 

*  * 

2 

■V 

** 

0 

!V 

** 

0 

!V 

** 

1 

!V 

*  * 

1 

!V 

** 

2 

!V 

*  * 

2 

■V 

** 

0 

!V 

*  * 

0 

!V 

** 

1 

!V 

** 

1 

!V 

** 

2 

!V 

*  * 

2 
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0514 
0515 
0516 
0517 
0518 
0519 
0520 
0521 
0522 
0523 
0524 
0525 
0526 
0527 
0528 
0529 
0530 
0531 
0532 
0533 
0534 
0535 
0536 
0537 
0538 
0539 
0540 
0541 
0542 
0543 
0544 
0545 
0546 
0547 
0548 
0549 
0550 
0551 
0552 
0553 
0554 
0555 
0556 
0557 
0558 
0559 
0560 
0561 
0562 
0563 
0564 
0565 
0566 
0567 
0568 
0569 
0570 


2020 

2021 
2022 
2023 
2025 
2026 
2499 


2498 

2024 

2500 
C 
C 
C 

C 
C 

C 


2036 


123 
2000 


2004 
2001 


ALPHA(1,4) 
LAMDA(1,4) 
ALPHA(2,4) 
LAMDA(2,4) 
ALPHA(3,4) 
LAMDA(3,4) 


4.5 

0. 

12. 

0. 

0. 

0. 


■V 

** 

0 

!V 

*  * 

0 

!V 

** 

1 

!V 

*  * 

1 

IV 

*  * 

2 

IV 

*  * 

2 

DO  2500  J  =  1,4 


WRITE 

WRITE 

FORMAT 

WRITE 

WRITE 

FORMAT 

WRITE 

FORMAT 

WRITE 

FORMAT 

WRITE 

FORMAT 

WRITE 

FORMAT 

WRITE 


,8A1) 


2,' 

) 

*    VP 

** 

0 

2,' 

) 

*    VP 

** 

1 

2,' 

) 

*    VP 

*• 

2 

2*) 


) 


5,2000) 

5,2020) (ITEXT(J,K) ,K=1,8) 

'0-2  CONSUMPTION  EQUATION  FOR 

5,2000) 

5,2021)ALPHA(1,J) ,LAMDA(1,J) 

'  VDOTO-2  =  ',F5.2,'  *  (  M  **  ',F5 

5,2022)ALPHA(2,J) ,LAMDA(2,J) 

+  ' ,F5.2, '  *  (  M  **  * ,F5 

5,2023)ALPHA(3,J) ,LAMDA(3,J) 

+  '  ,F5.2,*  *  (  M  **  •  ,F5 

5,2025) 

'  WHERE  VDOTO-2  IS  IN  ML( STP)/KG/MIN , 

5,2026) 

'  M  IS  IN  KG,  AND  VP  IS  IN  M/SEC  ' ) 

5,2000) 

CALL  PRMPTA( '  DO  YOU  WANT  TO  CHANGE  THIS  EXPRESSION? 
1        (Y  OR  N) ' ,IANS,NF) 
WRITE  (5,2000) 
IF  (NF.EQ.O)GO  TO  2498 
IF  (IANS(l) .EQ.'N* )GO  TO  2498 

IF  (IANS(l) .EQ. 'Y' )CALL  CHGEXP( ALPHA, LAMDA, ITEXT, J , LP) 
IF  (IANS(l) .EQ. ' Y' )GO  TO  2500 
GO  TO  2499 
WRITE(LP,2000) 

WRITE (LP, 2024) (ITEXT(J,K) ,K=1,8) 

FORMATdX,'  STANDARD  EXPRESSION  USED  FOR  0-2 
1  CONSUMPTION  FOR  ' , 8A1 ) 
CONTINUE 

MAIN  LOOP  FOR  ALL  MINERS 


DO  990  I  =  1,50 

STOP  SCREEN  EVERY  5  MINERS  TO  ALLOW  EXAMINATION  OF  THE  DATA 


IF  (MOD(I,5) .NE.0)GO  TO  123 

WRITE  (LP, 2000) 

WRITE  (5,2036) 

FORMAT ( '  "S"  TO  STOP,  RETURN  TO  CONTINUE') 

READ(5,2100)ANS 

IF  (ANS.EQ. 'S' )STOP 

WRITE  (LP, 2000) 

FORMAT ( '  ' ) 

IF  (JUNSTR(I) .EQ.0)GO  TO  1999 

WRITE  (LP, 2004)1, JUNSTR(I) 

FORMATdX,'  MINER  ',13,'  STARTING  AT  JUNCTION 

WRITE  (LP, 2001)1, MASS(I) ,VDOT(I) 

FORMATdX,'  MINER  ',13,'    MASS: 


F5.2 


KG 


13) 


VDOT  0-2 


,F5.2 
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0571  1        '  ML  0-2/KG/MIN') 

0572  ICURJUN  =  JUNSTR(I)  ! CURRENT  JUNCTION 
057  3             ICURPAT  =  ESCPAT(I)  ! CURRENT  PATH 

0574  IF(JUNSTR(I) .EQ.AIRJUN( ESCAIR( ESCPAT( I ) ,1) ,1) )GO  TO  137 

0575  IF(JUNSTR(I) . EQ.AIRJUN( ESCAIR( ESCPAT( I ) ,1) ,2) )GO  TO  137 

0576  WRITE  (LP, 2000) 

0577  WRITE  (LP, 2039) 

0578  WRITE  (LP, 2000) 

0579  2039     FORMAT ( '   !  I  !!!  I  !!!!!!!!!  I  !!!!  I  !!!!!!!!!!!  !  ' ) 

0580  WRITE  (LP,2032)ESCAIR(ESCPAT(I) ,1) ,JUNSTR(I) 

0581  2032     FORMAT ( '  MINER  IN  AIRWAY  ',13,'  THIS 

0582  1  AIRWAY  DOES  NOT  START  AT  JUNCTION  ',13) 

0583  WRITE  (5,2037) 

0584  2037     FORMAT( '  CONTINUE?  (  Y/N  ) ' ) 

0585  READ  (5,2100)ANS 

0586  2100     FORMAT(A) 

0587  IF  (ANS.EQ. 'N' ) STOP 

0588  C 

0589  C        NULL  0-2  USAGE 

0590  C 

0591  137      02USE  =  0. 

0592  C 

0593  C        NULL  TOTAL  TIME 

0594  C 

0595  TIMTOT  =  0. 

0596  C 

0597  C        NULL  ESCAPE  PATH  MARKER 

0598  C 

0599  INXT  =  0 

0600  C 

0601  C        NULL  REST  JUNCTION  FLAG  AND  COUNTER 

0602  C 

0603  IRESFLG  =  0 

0604  IRESCNT  =  0 

0605  C 

0606  C        NULL  TOTAL  DISTANCE  TRAVELED 

0607  C 

0608  SAIR=0. 

0609  C 

0610  1009     WRITE  ( LP, 2005 ) ICURPAT 

0611  2005     FORMAT  (IX,'  MINER  IS  ON  PATH  ',13) 

0612  DO  1000  IAIRJ  =  1,7 

0613  ICURAIR  =  ESCAIR( ICURPAT, IAIRJ) 

0614  IF  (ICURAIR. EQ.0)  GO  TO  989 

0615  IF    (IAIRJ. EQ.DGO   TO    1011 

0616  C 

0617  C        CHECK  TO  SEE  THAT  CURRENT  AIRWAY  AND  PREVIOUS  AIRWAY  INTERSECT 

0618  C 

0619  IF  (AIRJUN(ICURAIR,1) 

0620  1        .EQ. 

0621  2        AIRJUN(ESCAIR(ICURPAT,IAIRJ-1) ,2) )GO  TO  1011 

0622  IF  (AIRJUN(ICURAIR,1) 

0623  1        .EQ. 

0624  2        AIRJUN(ESCAIR( ICURPAT, IAIRJ-1) ,1) ) GO  TO  1011 

0625  IF  (AIRJUN(ICURAIR,2) 

0626  1        .EQ. 

0627  2        AIRJUN(ESCAIR( ICURPAT, IAIRJ-1) ,2) )GO  TO  1011 
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0628 

0629 

0630 

0631 

0632 

0633 

0634 

0635 

2033 

0636 

0637 

0638 

0639 

0640 

1011 

0641 

C 

0642 

C 

0643 

C 

0644 

0645 

0646 

0647 

C 

0648 

C 

0649 

C 

0650 

0651 

0652 

0653 

0654 

0655 

0656 

0657 

C 

0658 

C 

0659 

C 

0660 

C 

0661 

0662 

C 

0663 

C 

0664 

C 

0665 

0666 

C 

0667 

0668 

0669 

3006 

0670 

0671 

0672 

0673 

0674 

C 

0675 

C 

0676 

C 

0677 

141 

0678 

0679 

C 

0680 

0681 

c 

0682 

c 

0683 

c 

0684 

IF  (AIRJUN(ICURAIR,2) 

1  .EQ. 

2  AIRJUN(ESCAIR(ICURPAT,IAIRJ-1) ,1) )GO  TO  1011 
WRITE  (LP, 2000) 

WRITE  (LP, 2039) 

WRITE  (LP, 2000) 

WRITE  (LP,2033)ICURAIR,ESCAIR(ICURPAT,IAIRJ-1) 

FORMAT ( '  CURRENT  AIRWAY  #  *,I3,'  DOES  NOT  INTERSECT  WITH 

1   THE  PREVIOUS  AIRWAY  #  ',13) 

WRITE  (5,20  37) 

READ  (5,2100)ANS 

IF  (ANS.EQ.'N' )STOP 

WRITE  (LP, 2000) 

IS  THIS  A  REST  JUNCTION? 

IF  (AIRJUN(ESCAIR(ICURPAT,IAIRJ) ,1) . EQ.RESJUN( I , 1 ) )IRESFLG  =  1 
IF  (AIRJUN(ESCAIR(ICURPAT,IAIRJ) ,1) .EQ.RESJUN( I , 2 ) ) IRESFLG  =  2 
IF  (AIRJUN(ESCAIR(ICURPAT,IAIRJ) ,1) .EQ.RESJUN( I ,3 ) )IRESFLG  =  3 

REST  IF  APPROPRIATE 

IF(JFLG.EQ.l)  IRESFLG=0 

IF  ( IRESFLG. EQ.0) GO  TO  141 

IRESCNT  =  IRESCNT  +  1 

TIME  (21  +  IRESCNT)  =  RESTIM( I , IRESFLG) 

TIMTOT  =  TIMTOT  +  TIME(21  +  IRESFLG) 

W  =  MASS(I) 

VP  =  0.0  ! VELOCITY 

CALCULATE  0-2  UPTAKE 

USE  ONLY  THE  CONSTANT  TERM  FROM  THE  WALKING  EXPRESSION. 

U(21  +  IRESCNT)  =  ALPHA(1,3)  *  W  **  LAMDA(1,3) 

CALCULATE  0-2  USAGE 

02USE  =  02USE  +  U(21  +  IRESCNT)  *  W  *  TIME(21  +  IRESCNT)  /  1000 

WRITE  (5,2000) 

WRITE  (5,300  6)I,RESJUN(I,IRESFLG) ,RESTIM(I , IRESFLG) 

FORMATdX,'  MINER  ',13,'  IS  RESTING  AT  JUNCTION  ',13, 

1        *  FOR  '  ,F5.2,  *  MINUTES* ) 

WRITE  (5, 2002) TIMTOT 

WRITE  (5,2038)O2USE 

IRESFLG  =  0 

CALCULATE  VELOCITY  IN  AIRWAY  IAIRW 

H  =  AIRHGT(ICURAIR)      1AIRWAY  HEIGHT 
W  =  MASS(I)  ! MINER  MASS 

V=A+B*H+C*H**2+D*H**3 

CALCULATE  TIME  ON  PATH 

TIME(INXT  +  IAIRJ)  =  AIRLEN ( ICURAIR)/V 
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0685  TIMTOT  =  TIMTOT  +  TIME(INXT  +  IAIRJ) 

0686  C 

0687  C  CONVERT  VELOCITY  TO  M/SEC 

0688  C 

0689  C  FT      12  IN      0.0254  M      MIN 

0690  C  * * * 

0691  C  MIN      FT  IN         60  SEC 

0692  C 

0693  VP  =  V  *  12./60.*0.0254 

0694  C 

069  5  C  CALCULATE  0-2  UPTAKE 

0696  C 

0697  C  WE  HAVE  4  DIFFERENT  0-2  UPTAKE  EXPRESSIONS.   THEY  ARE: 

0698  C 

0699  C  CRAWL  H  <  30" 

0700  C  VDOT02  =  48.6  *  VP 

0701  C  LOCMTN  =  1 

0702  C 

0703  C  DUCKWALK         30"  <  H  <  50" 

0704  C  VDOT02  =  39.6  *  VP 

0705  C  LOCMTN  =  2 

0706  C 

0707  C  WALK  50"  <  H 

0708  C  VDOT02  =  7.06  +3.6  *  VP  *  VP 

0709  C  LOCMTN  =  3 

0710  C 

0711  C  RUN  S  >  134  M/MIN 

0712  C  VDOT02  =  4.5  +  12.VP 

0713  C  LOCMTN  =  4 

0714  C 

0715  IF  (H.LE.30.0)LOCMTN  =  1 

0716  IF  (H.GT.30.0.AND.H.LE.50.0)LOCMTN  =  2 

0717  IF  (H.GT.50.0.AND.H.LE.84.0)LOCMTN  =  3 

0718  IF  (H.GT. 84.0. AND. IWORR.EQ. 1 ) LOCMTN  =  4 

0719  IF  ( LOCMTN. EQ. 4 )VP  =  V5MPH 

0720  U(INXT  +  IAIRJ)  =  0. 

0721  DO  149  J  =  1,3 

07  22  U(INXT  +  IAIRJ)  =  U(INXT  +  IAIRJ)  + 

0723  1        ALPHA( J, LOCMTN)  *  W  **  LAMDA( J , LOCMTN)  *  VP  **  (J-l) 

0724  149  CONTINUE 

0725  C 

0726  C  CALCULATE  0-2  USAGE 

0727  C 

0728  153  02USE  =  02USE  +  U(INXT  +  IAIRJ)  *  W  *  TIME(INXT  +  IAIRJ)  /  1000 

0729  C 

0730  WRITE  (LP, 2000) 

0731  WRITE  (LP, 2006)1, ICURAIR 

0732  2006     FORMATQX,'  MINER  ',13,'  IS  IN  AIRWAY  ',13) 

0733  WRITE  ( LP, 2027 )H , ( ITEXT( LOCMTN , J) , J=l ,8 ) 

0734  2027     FORMATdX,'  AIRWAY  HEIGHT  =  ',F5.2,'  IN.   LOCOMOTION  =>  ',8A1) 

0735  WRITE  (LP, 2002)  TIMTOT 

0736  2002     FORMATdX,'  TIME  USED  SINCE  START  OF  ESCAPE :', F6 . 3 , '  MINUTES') 

0737  WRITE  ( LP, 2038 )02USE 

0738  2038     FORMATdX,'  0-2  USED  SINCE  START  OF  ESCAPE:  ',F8.2,'  1  STP') 

0739  IF(TIMTOT.GT. 60.0. OR. 02USE.GT. 100.0)  WRITE( LP, 7700 ) 

0740  7700     FORMAT(/5X,'  *****    MINER  STOPS   *****•/) 

0741  IJI=INXT+IAIRJ 
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0742 

0743 

2041 

0744 

0745 

0746 

0747 

0748 

2042 

0749 

0750 

0751 

0752 

0753 

2046- 

0754 

0755 

C 

0756 

0757 

0758 

0759 

C 

0760 

C 

0761 

C 

0762 

0763 

0764 

0765 

0766 

0767 

0768 

0769 

0770 

0771 

0772 

0773 

0774 

0775 

0776 

0777 

0778 

2034 

0779 

0780 

0781 

0782 

0783 

1111 

0784 

0785 

1000 

0786 

c 

0787 

c 

0788 

c 

0789 

c 

0790 

c 

0791 

989 

0792 

C 

0793 

C 

0794 

C 

0795 

0796 

140 

0797 

C 

0798 

C 

WRITE(LP,2041)  I , ICURAIR,U( IJI ) 

FORMATQX,'  0-2  UPTAKE  RATE  BY  MINER'  ,  IX,  13  ,  IX,  '  IN  AIRWAY', 

1   IX, 13, IX, 'IS  ' ,1X,F5.2,1X,'ML  0-2  /KG  /MIN' ) 

VS=V/88. 

SAIR=SAIR+AIRLEN( ICURAIR) 

WRITE(LP,2042)  V,VS ,AIRLEN( ICURAIR) ,SAIR 

F0RMAT(1X,'  MINER  SPEED* , 1X,F7 . 2 , IX, ' FPM ' ,4X,' (  ',F7.2,'  MPH) * 

1  ,4X,' AIRWAY  LENGTH', 

2  IX, F7. 2, IX, 'FT* ,4X, 

3  /,2X, 'TOTAL  DISTANCE  TRAVELLED: ', IX, F7 . 2 , IX, ' FT ' ) 
WRITE (LP, 2046)  I,VDOT(I) 

FORMATdX,*  MAX  0-2  UPTAKE  BY  MINER  '  ,  IX,  13  ,  IX,  '  IS  ',F5.2 
1  ,1X'ML  0-2/KG/MIN') 

IF  (MOD(IAIRJ,7) .NE.0)GO  TO  1000 

IF  (NXTPAT(ICURPAT) .EQ.0)GO  TO  1000 

INXT  =  INXT  +  7  ! MARKER  FOR  NUMBER  OF  ESCAPE  PATHS 

CHECK  TO  SEE  THAT  CURRENT  AIRWAY  AND  NEXT  AIRWAY  INTERSECT 

IF  (AIRJUN( ICURAIR, 2) 

1  .EQ. 

2  AIRJUN(ESCAIR(NXTPAT(ICURPAT) ,1) ,1) ) GO  TO  1111 
IF  (AIRJUN( ICURAIR, 2) 

1  .EQ. 

2  AIRJUN(ESCAIR(NXTPAT(ICURPAT) ,1) ,2) ) GO  TO  1111 
IF  (AIRJUN( ICURAIR, 1) 

1  .EQ. 

2  AIRJUN(ESCAIR(NXTPAT(ICURPAT) ,1) ,1) )GO  TO  1111 
IF  (AIRJUN( ICURAIR, 1) 

1  .EQ. 

2  AIRJUN(ESCAIR(NXTPAT(ICURPAT) ,1) ,2) )GO  TO  1111 
WRITE  (LP, 2000) 

WRITE  (LP, 2039) 

WRITE  (LP, 2000) 

WRITE  (LP,2034)ICURAIR,ESCAIR(NXTPAT(ICURPAT) ,1) 

FORMAT ( '  CURRENT  AIRWAY  #  ',13,'  DOES  NOT  INTERSECT  WITH 

1   THE  NEXT  AIRWAY  #  ',13) 

WRITE  (5,2037) 

READ(5,2100)ANS 

IF  (ANS.EQ. 'N' )STOP 

ICURPAT  -  NXTPAT(ICURPAT) 

GO  TO  1009 

CONTINUE 

CALCULATE  PHYSIOLOGICAL  STRESS  OF  THIS  WORK 

NULL  G  FACTOR 

G  =  0. 

CALCULATE  G  FACTOR 

DO  140  ITOT  =  1,INXT  +  IAIRJ  -1 

G  =  G  +  U(ITOT)  *  TIME(ITOT)  /  VDOT ( I )  /  TIMTOT 

ADD  FOR  REST  TIMES 
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0799  C 

0800  DO  142  IRES  =  1,IRESCNT 

0801  142  G  =  G  +  U(IRES)  *  TIME(21  +  IRES)  /  VDOT(I)  /  TIMTOT 

0802  C 

0803  C  CALCULATE  TAU2 

0804  C  THIS  IS  TIME  AVAILABLE  TO  A  MINER  AT  MAXIMUM  PHYSICAL 

0805  C  EXERTION  BEFORE  MINER  STOPS 

0806  C 

0807  TAU2  =  10.  **  (3.76  -  3.1  *  G) 

0808  WRITE  (LP,2003)TAU2 

0809  2003  FORMATdX,'  TIME  AVAILABLE  FOR  ESCAPE:  ',G10.2,'  MINUTES  ') 

0810  990  CONTINUE 

0811  1999  CLOSE  (UNIT  =  FDAT) 

0812  STOP 

0813  END 
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0001 

C 

0002 

C 

0003 

C 

0004 

C 

0005 

C 

0006 

C 

0007 

C 

0008 

C 

0009 

C 

0010 

C 

0011 

C 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 

0021 

0022 

0023 

ROUTINE  TO  PUT  PROMPT  STRING  OUT  TO  TERMINAL 
AND  THEN  GET  BACK  FROM  TERMINAL  A  STRING  AND 
ITS  LENGTH.  THIS  IS  VERY  USEFUL  FOR  GETTING 
FILE  NAMES  FROM  THE  TERMINAL. 

IN   >>  INPUT  STRING  FROM  CALLING  ROUTINE 
OUT  >>  STRING  FROM  TERMINAL  BACK  TO 

CALLING  ROUTINE. 
M    >>  LENGTH  OF  OUT. 

SUBROUTINE  PRMPTA( IN ,OUT ,M) 

BYTE  IN(80) ,OUT(80) 

DO  10  1=1,80 

IF(IN(I)  ,EQ.  0)GOTO  11 

10  CONTINUE 

11  K=I-1 

WRITE (5, 1000) (IN(I) ,I=1,K) 
1000    FORMAT ( IX, <K>A1,$) 

READ(5,2000)  M, (OUT( L) , L=l ,M) 
2000    FORMAT ( Q, 8 0A1) 

RETURN 

END 
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0001  C 

0002  C  ROUTINE  TO  PUT  PROMPT  STRING  OUT  TO  TERMINAL 

000  3  C        AND  THEN  GET  BACK  FROM  TERMINAL  A  FLOATING  POINT  VALUE 

0004  C        AND  ITS  LENGTH.   THIS  IS  VERY  USEFUL  FOR  GETTING 

000  5  C        VALUES  FROM  THE  TERMINAL  WHERE  A  CARRIAGE 

0006  C  RETURN  MAY  BE  EXPECTED. 

0007  C 

0008  C  IN   >>  INPUT  STRING  FROM  CALLING  ROUTINE 

0009  C  XID  >>  VALUE  FROM  TERMINAL  TRANSFERED  BACK 

0010  C  TO  CALLING  ROUTINE.   XID  IS  READ  AS  A 

0011  C  STRING  AND  DECODED  BEFORE  BEING  SENT  BACK 

0012  C  NZ   >>  LENGTH  OF  OUT,  NZ  =  0  IF  CARRIAGE 

0013  C  RETURN  IS  THE  RESPONSE  FROM  THE  TERMINAL. 

0014  C 

0015  SUBROUTINE  PRMPTF( IN,XID,NZ ) 

0016  BYTE  IN(80) ,BUFF(80) 

0017  CALL  ERRSET( 64, .TRUE. , .FALSE. , .TRUE. , .FALSE. ,200) 

0018  DO  10  1=1,80 

0019  IF(IN(I)  .EQ.  0)GOTO  11 

0020  10  CONTINUE 

0021  11  K=I-1 

0022  12  WRITE(5,900) 

0023  900     FORMAT(IX) 

0024  WRITE(5,1000) (IN(L) ,L=1,K) 

0025  1000    F0RMAT(1X,<K>A1,1X,$) 

0026  READ(5,2002)NZ,(BUFF(KL) ,KL=1,NZ) 

0027  2002    FORMAT(Q,20A1) 

0028  DO  40  M=1,NZ 

0029  IF(BUFF(M)  .NE.  ',')GOTO  40 

0030  NY=M 

0031  GOTO  30 

0032  40  CONTINUE 

0033  DECODE (NZ, 2000, BUFF, ERR=13)X 

0034  2000    FORMAT(F8.0) 

0035  XID=X 

0036  RETURN 

0037  13  CONTINUE 

0038  DO  20  1=1, NZ 

0039  IF(BUFF(I)  .EQ.  '.')  GOTO  20 

0040  IF(BUFF(I)  .EQ.  '  ')  GOTO  20 

0041  IF(BUFF(I)  .GE.  "60  .AND.  BUFF(I)  .LE.  "71)  GOTO  20 

0042  IF(BUFF(I)  .EQ.  0 ) STOP 

0043  NY=I 

0044  GOTO  30 

0045  20  CONTINUE 

0046  30  MX=NY+K+1 

0047  WRITE(5,1010) 

0048  1010    FORMAT(<MX>X, ' ~ • ) 

0049  WRITE(5,1001) 

0050  1001    FORMAT (IX,'  ERROR   RE-ENTER  DATA') 

0051  GOTO  12 

0052  END 
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0001 

C 

0002 

C 

0003 

C 

0004 

C 

0005 

0006 

C 

0007 

0008 

c 

0009 

0010 

c 

0011 

2000 

0012 

1 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

2 

0020 

0021 

100 

0022 

0023 

0024 

0025 

0026 

0027 

0028 

10 

0029 

0030 

3 

0031 

0032 

200 

0033 

0034 

0035 

0036 

0037 

0038 

20 

0039 

0040 

21 

0041 

0042 

2020 

0043 

0044 

0045 

2021 

0046 

0047 

2022 

0048 

0049 

2023 

0050 

0051 

2025 

0052 

0053 

2026 

0054 

0055 

0056 

SUBROUTINE  TO  CHANGE  COEFFICIENT  VALUES  FOR  0-2 
CONSUMPTION  EQUATIONS. 

SUBROUTINE  CHGEXP( ALPHA, LAMDA, ITEXT, J, LP) 

REAL  *4  ALPHA(3,4) ,LAMDA(3,4) 

BYTE  IANS(5) , ITEXT (4, 8) 

FORMAT ('  ') 

CALL  PRMPTA( '  CHANGE  COEFFICIENT  OR  EXPONENT?  (C  OR  E)  RETURN 

1  TO  STOP  ' ,IANS,NF) 

WRITE  (5,2000) 

IF  (NF.EQ.O)GO  TO  21 

IF  (IANS(l) .EQ.'C )GO  TO  2 

IF  (IANS(l) .EQ.'E' )GO  TO  3 

GO  TO  1 

DO  10  I  =  1,3 

WRITE (5, 100) (1-1) ,ALPHA(I,J) 

FORMAT ( '  COEFFICIENT  FOR  VP  **  ',12,'  =  ',F5.2) 

CALL  PRMPTF( '  NEW  VALUE?  RETURN  FOR  NO  CHANGE  *,VNEW,NF) 

WRITE  (5,2000) 

IF  (NF.NE.0)ALPHA(I,J)  =  VNEW 

WRITE  (5,2000) 

WRITE  (5,100)(I-1),ALPHA(I,J) 

WRITE  (5,2000) 

CONTINUE 

GO  TO  1 

DO  20  I  =  1,3 

WRITE(5,200)(I-1) ,LAMDA(I,J) 

FORMAT('  EXPONENT  FOR  VP  **  ',12,'  =  ',F5.2) 

CALL  PRMPTF('  NEW  VALUE?  RETURN  FOR  NO  CHANGE  ',VNEW,NF) 

IF  (NF.NE.0)LAMDA(I,J)  =  VNEW 

WRITE  (5,2000) 

WRITE  (5,200) (1-1) ,LAMDA(I,J) 

WRITE  (5,2000) 

CONTINUE 

GO  TO  1 

WRITE(LP,2000) 

WRITE (LP, 20 20) ( ITEXT ( J,K) ,K=1 ,8 ) 

FORMAT ( '  0-2  CONSUMPTION  EQUATION  FOR  ',8A1) 

WRITE(LP,2000) 

WRITE ( LP , 20  2 1 ) ALPHA( 1 , J ) , LAMDA( 1 , J ) 

FORMAT ( '    VDOTO-2    =    ',F5.2,'     *     (    M    **     ',F5.2,'     )     *    VP    **    0') 

WRITE(LP,2022)ALPHA(2,J) ,LAMDA(2,J) 

FORMATC  +    ',F5.2,'     *     (    M    **     ',F5.2,'     )     *    VP    **    1') 

WRITE(LP,20  23)ALPHA(3,J) ,LAMDA(3,J) 

FORMATC  +    ',F5.2,'     *     (    M    **     ',F5.2,'     )     *    VP    **    2') 

WRITE(LP,2025) 

FORMAT ( '    WHERE    VDOTO-2    IS    IN    ML( STP)/KG/MIN , ' ) 

WRITE(LP,2026) 

FORMATC  M  IS  IN  KG,  AND  VP  IS  IN  M/SEC  ') 

WRITE(LP,2000) 

RETURN 

END 
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APPENDIX  D 

.— LISTIN 

G  OF  INI 

1 

1 

6 

20. 

600. 

3 

2 

7 

25. 

600. 

5 

3 

8 

30. 

600. 

7 

4 

9 

35. 

600. 

9 

5 

10 

40. 

600. 

11 

6 

11 

45. 

600. 

13 

7 

12 

50. 

600. 

15 

8 

13 

55. 

600. 

17 

9 

14 

60. 

600. 

19 

10 

15 

65. 

600. 

21 

11 

16 

70. 

600. 

23 

12 

17 

75. 

600. 

25 

13 

18 

80. 

600. 

27 

14 

19 

77. 

600. 

29 

15 

20 

72. 

600. 

31 

16 

21 

67. 

600. 

33 

17 

22 

62. 

600. 

35 

18 

23 

57. 

600. 

37 

19 

24 

52. 

600. 

39 

20 

25 

47. 

600. 

2 

1 

2 

22. 

600. 

4 

2 

3 

27. 

600. 

6 

3 

4 

32. 

600. 

8 

4 

5 

37. 

600. 

10 

6 

7 

42. 

600. 

12 

7 

8 

47. 

600. 

14 

8 

9 

52. 

600. 

16 

9 

10 

57. 

600. 

18 

11 

12 

62. 

600. 

20 

12 

13 

67. 

600. 

22 

13 

14 

72. 

600. 

24 

14 

15 

77. 

600. 

26 

16 

17 

80. 

600. 

28 

17 

18 

75. 

600. 

30 

18 

19 

70. 

600. 

32 

19 

20 

65. 

600. 

34 

21 

22 

60. 

600. 

36 

22 

23 

55. 

600. 

38 

23 

24 

50. 

600. 

40 

24 

25 

45. 

600. 

41 

13 

26 

80. 

600. 

999 

999 

999 

999 

0. 

THIS  IS  AN  AIRWAY  LIST  FOR  MINE  1.   IT  GETS  TRANSLATED  BY  THE 
PROGRAM  TRANS  INTO  A  BINARY  FILE.   THE  FORMAT  IS: 

AIRWAY   JUN      JUN      HEIGHT   LENGTH 
NUMBER   1        2        (IN.)    (FT.) 

JUNCTION  1  IS  THE  STARTING  JUNCTION  FOR  THE  AIRWAY  AND 
JUNCTION  2  IS  THE  ENDING  JUNCTION  FOR  THE  AIRWAY. 

ENTER  0  FOR  REMAINING  JUNCTIONS  IF  THERE  ARE  LESS  THAN  6  JUNCTIONS 
IN  AN  AIRWAY. 

AIRWAY  #  999  STOPS  THE  PROGRAM. 
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1 

1 

1 

34.1 

1 

78.6 

3 

5.       13 

5. 

0 

0. 

2 

2 

2 

44.52634 

2 

86. 

3 

5.       13 

5. 

0 

0. 

3 

3 

3 

24.51376 

3 

87. 

13 

5.       0 

0. 

0 

0. 

4 

4 

4 

28.21625 

4 

88. 

3 

5.       13 

5. 

0 

0. 

5 

5 

5 

35.38699 

5 

89. 

3 

5.       13 

5. 

0 

0. 

6 

6 

6 

42.54644 

6 

90. 

8 

5.       13 

5. 

0 

0. 

7 

7 

7 

22.44606 

7 

91. 

8 

5.       13 

5. 

0 

0. 

8 

8 

8 

45.57353 

8 

92. 

13 

5.       0 

0. 

0 

0. 

9 

9 

9 

40.51924 

9 

93. 

8 

5.       13 

5. 

0 

0. 

10 

10 

10 

42.15463 

10 

94. 

8 

5.       13 

5. 

0 

0. 

11 

11 

11 

42.9893 

11 

95. 

13 

5.       0 

0. 

0 

0. 

12 

12 

12 

28.65271 

12 

53. 

13 

5.       0 

0. 

0 

0. 

13 

13 

13 

35.98224 

13 

84. 

0 

0.       0 

0. 

0 

0. 

14 

14 

14 

27.58545 

14 

83. 

13 

5.       0 

0. 

0 

0. 

15 

15 

15 

34.76538 

15 

82. 

13 

5.       0 

0. 

0 

0. 

16 

16 

16 

40.26243 

16 

81. 

18 

300      13 

5. 

0 

0. 

17 

17 

17 

42.99278 

17 

80. 

18 

5.       13 

5. 

0 

0. 

18 

18 

18 

31.3226 

18 

79. 

13 

5.       0 

0. 

0 

0. 

19 

19 

19 

32.55729 

19 

78. 

18 

5.       13 

5. 

0 

0. 

20 

20 

20 

33.4361 

20 

77. 

18 

5.       13 

5. 

0 

0. 

21 

21 

21 

40.95516 

21 

76. 

23 

5.       13 

5. 

0 

0. 

22 

22 

22 

37.58183 

22 

75. 

23 

5.       13 

5. 

0 

0. 

23 

23 

23 

22.83065 

23 

74. 

13 

5.       0 

0. 

0 

0. 

24 

24 

24 

44.7841 

24 

73. 

23 

5.       13 

5. 

0 

0. 

25 

25 

25 

22.73881 

25 

72. 

23 

5.       13 

5. 

0 

0. 

26 

1 

1 

22.5 

26 

71.2 

3 

5.       13 

5. 

0 

0. 

27 

5 

5 

30.8 

27 

77.6 

3 

5.       13 

5. 

0 

0. 

28 

6 

6 

36.8 

28 

78. 

8 

5.       13 

5. 

0 

0. 

29 

10 

10 

33.5 

38 


29 

80.7 

8 

5. 

13 

30 

17 

17 

38.9 

30 

81.3 

18 

5. 

13 

31 

20 

20 

42. 

31 

82.9 

18 

5. 

13 

32 

25 

25 

34.1 

32 

78.6 

23 

5. 

13 

999 

999 

999 

999 

THIS  IS  A  MINER  LIST  FOR  MINE  1.   IT  GETS  TRANSLATED  BY  THE 
PROGRAM  TRANS  INTO  A  BINARY  FILE.   THE  FORMAT  FOR  THE  FIRST  LINE  IS 


MINER 
NUMBER 


START 
JUNCTION 


ESCAPE 
PATH 


MAX  0-2 

INTAKE 


MINER  #  999  STOPS  THE  PROGRAM. 


THIS  IS  A  MINER  MASS  LIST  FOR  MINE  1.   IT  GETS  TRANSLATED  BY  THE 
PROGRAM  TRANS  INTO  A  BINARY  FILE.   THE  FORMAT  FOR  THE  SECOND  LINE  IS: 


MINER 

MINER 

REST 

REST 

REST 

REST 

REST 

NUMBER 

MASS 

JUNCT 

TIME 

JUNCT 

TIME 

TIME 

1 

1 

2 

2 

3 

MINER  #  999  STOPS  THE  PROGRAM. 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

26 

21 

22 

23 

24 

25 

999 


2 

4 

5 

15 

4 

5 

15 

41 

5 

15 

41 

0 

6 

5 

15 

41 

8 

6 

5 

15 

10 

12 

15 

41 

12 

15 

41 

0 

15 

41 

0 

0 

14 

15 

41 

0 

16 

14 

15 

41 

18 

20 

41 

0 

20 

41 

0 

0 

41 

0 

0 

0 

22 

41 

0 

0 

24 

22 

41 

0 

26 

28 

25 

41 

28 

25 

41 

0 

25 

41 

0 

0 

30 

25 

41 

0 

32 

27 

17 

7 

1 

11 

18 

20 

34 

36 

35 

25 

36 

35 

25 

41 

35 

25 

41 

0 

38 

35 

25 

41 

40 

38 

35 

25 

999 

999 

999 

0 

41 

0 

0 

0 

41 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

6 

41 

41 

0 

0 

0 

41 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

26 

0 

0 

0 

0 

0 

0 

0 


THIS  IS  A  ESCAPE  PATH  LIST  FOR  MINE  1.   IT  GETS  TRANSLATED  BY  THE 
PROGRAM  TRANS  INTO  A  BINARY  FILE.   THE  FORMAT  IS: 

PATH     AIRWAY   AIRWAY   AIRWAY   AIRWAY   AIRWAY   AIRWAY   AIRWAY   NEXT 
NUMBER   12        3        4        5        6        7        PATH 

IF  THE  PATH  REACHES  THE  EXIT  BEFORE  AIRWAY  7,  THEN  ENTER  0  FOR 

REMAINING  AIRWAYS.   IF  AIRWAY  7  DOES  NOT  TAKE  THE  PATH  TO  THE 

EXIT,  THE  PATH  MAY  BE  CONTINUED  BY  THE  PATH  NAMED  IN  THE  LAST  COLUMN. 


PATH  #  9  99  STOPS  THE  PROGRAM. 


40 


APPENDIX  E. — LISTING  OF  OUTPUT  DATA 

CURRENT  VALUE  OF  JFLG  IS      1 
CURRENT  VALUE  OF  IWORR  IS      0 

A  =  252.2 
B  =  -16.14 
C  =  0.4080 
D  =  -.2440E-02 

0-2  CONSUMPTION  EQUATION  FOR  CRAWLING 

VDOTO-2  =  0.00  *  (  M  **  0.00  )  *  VP  **  0 
+  48.60  *  (  M  **  0.00  )  *  VP  **  1 
+   0.00  *  (  M  **   0.00  )  *  VP  **  2 

WHERE  VDOTO-2  IS  IN  ML( STP)/KG/MIN , 

M  IS  IN  KG,  AND  VP  IS  IN  M/SEC 

0-2  CONSUMPTION  EQUATION  FOR  DUCKWALK 

VDOTO-2  =  0.00  *  (  M  **  0.00  )  *  VP  **  0 
+  39.60  *  (  M  **  0.00  )  *  VP  **  1 
+   0.00  *  (  M  **   0.00  )  *  VP  **  2 

WHERE  VDOTO-2  IS  IN  ML( STP)/KG/MIN , 

M  IS  IN  KG,  AND  VP  IS  IN  M/SEC 

0-2  CONSUMPTION  EQUATION  FOR  WALKING 

VDOTO-2  =   7.06  *  (  M 
+   0.00  *  (  M 
+   3.60  *  (  M 
WHERE  VDOTO-2  IS  IN  ML( STP ) /KG/MIN , 
M  IS  IN  KG,  AND  VP  IS  IN  M/SEC 

0-2  CONSUMPTION  EQUATION  FOR  RUNNING 

VDOTO-2  =   4.50  *  (  M 

+  12.00  *  (  M 

+   0.00  *  (  M 
WHERE  VDOTO-2  IS  IN  ML( STP)/KG/MIN , 
M  IS  IN  KG,  AND  VP  IS  IN  M/SEC 


MINER    1  STARTING  AT  JUNCTION    1 

MINER    1    MASS:  78.60  KG    VDOT  0-2:  34.10  ML  O-2/KG/MIN 

MINER  IS  ON  PATH    1 

MINER    1  IS  IN  AIRWAY    2 

AIRWAY  HEIGHT  =  22.00  IN.   LOCOMOTION  =>  CRAWLING 

TIME  USED  SINCE  START  OF  ESCAPE:   8.745  MINUTES 

0-2  USED  SINCE  START  OF  ESCAPE:     11.64  1  STP 

0-2  UPTAKE  RATE  BY  MINER    1  IN  AIRWAY    2  IS   16.94  ML  0-2  /KG  /MIN 

MINER  SPEED    68.61  FPM     (     0.78  MPH)     AIRWAY  LENGTH   600.00  FT 

TOTAL  DISTANCE  TRAVELLED:   600.00  FT 

MAX  0-2  UPTAKE  BY  MINER     1  IS  34.10  ML  O-2/KG/MIN 


** 

0.00 

) 

* 

VP 

**    0 

*  * 

0.00 

) 

* 

VP 

**    I 

*  * 

0.00 

) 

* 

VP 

**    2 

** 

0.00 

) 

*    VP 

**    0 

*  * 

0.00 

) 

*    VP 

**    2 

** 

0.00 

) 

*    VP 

**    2 
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MINER    1  IS  IN  AIRWAY    4 

AIRWAY  HEIGHT  =  27.00  IN.   LOCOMOTION  =>  CRAWLING 

TIME  USED  SINCE  START  OF  ESCAPE:  17.859  MINUTES 

0-2  USED  SINCE  START  OF  ESCAPE:     23.29  1  STP 

0-2  UPTAKE  RATE  BY  MINER    1  IN  AIRWAY    4  IS   16.25  ML  0-2  /KG  /MIN 

MINER  SPEED    65.83  FPM     (     0.75  MPH)     AIRWAY  LENGTH   600.00  FT 

TOTAL  DISTANCE  TRAVELLED:  1200.00  FT 

MAX  0-2  UPTAKE  BY  MINER     1  IS  34.10  ML  O-2/KG/MIN 

MINER    1  IS  IN  AIRWAY    5 

AIRWAY  HEIGHT  =  30.00  IN.   LOCOMOTION  =>  CRAWLING 

TIME  USED  SINCE  START  OF  ESCAPE:  26.513  MINUTES 

0-2  USED  SINCE  START  OF  ESCAPE:     34.93  1  STP 

0-2  UPTAKE  RATE  BY  MINER    1  IN  AIRWAY    5  IS   17.12  ML  0-2  /KG  /MIN 

MINER  SPEED    69.33  FPM     (     0.79  MPH)     AIRWAY  LENGTH   600.00  FT 

TOTAL  DISTANCE  TRAVELLED:  1800.00  FT 

MAX  0-2  UPTAKE  BY  MINER     1  IS  34.10  ML  O-2/KG/MIN 

MINER    1  IS  IN  AIRWAY   15 

AIRWAY  HEIGHT  =  55.00  IN.   LOCOMOTION  =>  WALKING 

TIME  USED  SINCE  START  OF  ESCAPE:  29.625  MINUTES 

0-2  USED  SINCE  START  OF  ESCAPE:     37.50  1  STP 

0-2  UPTAKE  RATE  BY  MINER    1  IN  AIRWAY   15  IS   10.51  ML  0-2  /KG  /MIN 

MINER  SPEED   192.78  FPM     (     2.19  MPH)     AIRWAY  LENGTH   600.00  FT 

TOTAL  DISTANCE  TRAVELLED:  2400.00  FT 

MAX  0-2  UPTAKE  BY  MINER     1  IS  34.10  ML  O-2/KG/MIN 

MINER    1  IS  IN  AIRWAY   41 

AIRWAY  HEIGHT  =  80.00  IN.   LOCOMOTION  =>  WALKING 

TIME  USED  SINCE  START  OF  ESCAPE:  31.483  MINUTES 

0-2  USED  SINCE  START  OF  ESCAPE:     39.95  1  STP 

0-2  UPTAKE  RATE  BY  MINER    1  IN  AIRWAY   41  IS   16.75  ML  0-2  /KG  /MIN 

MINER  SPEED   322.98  FPM     (     3.67  MPH)     AIRWAY  LENGTH   600.00  FT 

TOTAL  DISTANCE  TRAVELLED:  3000.00  FT 

MAX  0-2  UPTAKE  BY  MINER     1  IS  34.10  ML  O-2/KG/MIN 

TIME  AVAILABLE  FOR  ESCAPE:    0.20E+03  MINUTES 

MINER   32  STARTING  AT  JUNCTION   25 

MINER   32    MASS:  78.60  KG    VDOT  0-2:  34.10  ML  O-2/KG/MIN 

MINER  IS  ON  PATH   25 


MINER   32  IS  IN  AIRWAY   40 

AIRWAY  HEIGHT  =  45.00  IN.   LOCOMOTION  =>  DUCKWALK 

TIME  USED  SINCE  START  OF  ESCAPE:   4.623  MINUTES 

0-2  USED  SINCE  START  OF  ESCAPE:      9.49  1  STP 

0-2  UPTAKE  RATE  BY  MINER   32  IN  AIRWAY   40  IS   26.11  ML  0-2  /KG  /MIN 

MINER  SPEED   129.78  FPM     (     1.47  MPH)     AIRWAY  LENGTH   600.00  FT 

TOTAL  DISTANCE  TRAVELLED:   600.00  FT 

MAX  0-2  UPTAKE  BY  MINER    32  IS  34.10  ML  O-2/KG/MIN 

MINER   32  IS  IN  AIRWAY   38 

AIRWAY  HEIGHT  =  50.00  IN.   LOCOMOTION  =>  DUCKWALK 
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TIME  USED  SINCE  START  OF  ESCAPE:   8.368  MINUTES 

0-2  USED  SINCE  START  OF  ESCAPE:     18.97  1  STP 

0-2  UPTAKE  RATE  BY  MINER   32  IN  AIRWAY   38  IS   32.23  ML  0-2  /KG  /MIN 

MINER  SPEED   160.23  FPM     (     1.82  MPH)     AIRWAY  LENGTH   600.00  FT 

TOTAL  DISTANCE  TRAVELLED:  1200.00  FT 

MAX  0-2  UPTAKE  BY  MINER    32  IS  34.10  ML  0-2/KG/MIN 

MINER   32  IS  IN  AIRWAY   35 

AIRWAY  HEIGHT  =  57.00  IN.   LOCOMOTION  =>  WALKING 

TIME  USED  SINCE  START  OF  ESCAPE:  11.281  MINUTES 

0-2  USED  SINCE  START  OF  ESCAPE:     21.49  1  STP 

0-2  UPTAKE  RATE  BY  MINER   32  IN  AIRWAY   35  IS   11.00  ML  0-2  /KG  /MIN 

MINER  SPEED   205.98  FPM     (     2.34  MPH)     AIRWAY  LENGTH   600.00  FT 

TOTAL  DISTANCE  TRAVELLED:  1800.00  FT 

MAX  0-2  UPTAKE  BY  MINER    32  IS  34.10  ML  O-2/KG/MIN 

MINER   32  IS  IN  AIRWAY   25 

AIRWAY  HEIGHT  =  80.00  IN.   LOCOMOTION  =>  WALKING 

TIME  USED  SINCE  START  OF  ESCAPE:  13.138  MINUTES 

0-2  USED  SINCE  START  OF  ESCAPE:     23.94  1  STP 

0-2  UPTAKE  RATE  BY  MINER   32  IN  AIRWAY   25  IS   16.75  ML  0-2  /KG  /MIN 

MINER  SPEED   322.98  FPM     (     3.67  MPH)     AIRWAY  LENGTH   600.00  FT 

TOTAL  DISTANCE  TRAVELLED:  2400.00  FT 

MAX  0-2  UPTAKE  BY  MINER    32  IS  34.10  ML  O-2/KG/MIN 

MINER   32  IS  IN  AIRWAY   41 

AIRWAY  HEIGHT  =  80.00  IN.   LOCOMOTION  =>  WALKING 

TIME  USED  SINCE  START  OF  ESCAPE:  14.996  MINUTES 

0-2  USED  SINCE  START  OF  ESCAPE:     26.38  1  STP 

0-2  UPTAKE  RATE  BY  MINER   32  IN  AIRWAY   41  IS   16.75  ML  0-2  /KG  /MIN 

MINER  SPEED   322.98  FPM     (     3.67  MPH)     AIRWAY  LENGTH   600.00  FT 

TOTAL  DISTANCE  TRAVELLED:  3000.00  FT 

MAX  0-2  UPTAKE  BY  MINER    32  IS  34.10  ML  O-2/KG/MIN 

TIME  AVAILABLE  FOR  ESCAPE:     53.      MINUTES 


<rU.S.  CPO:     1985-605-017/20,146  int.-bu.of  mines,pgh.,p a.  28  194 
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U.S.  Department  of  the  Interior 
Bureau  of  Mines— Prod,  and  Distr. 
Cochrans  Mill  Road 
P.O.  Box  1807C 
Pittsburgh.  Pa.  15236 


OFFICIAL  BUSINESS 
PENALTY  FOR  PRIVATE  USE,  MOO 


AN  EQUAL  OPPORTUNITY  EMPLOYER 


POSTAGE  ANO  FEES  PAID 
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